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Abstract. This article reviews the recent theoretical and experimental advances 
in the study of ultracold gases made of bosonic particles interacting via the long- 
range, anisotropic dipole-dipole interaction, in addition to the short-range and 
isotropic contact interaction usually at work in ultracold gases. The specific 
properties emerging from the dipolar interaction are emphasized, from the mean- 
field regime valid for dilute Bose-Einstein condensates, to the strongly correlated 
regimes reached for dipolar bosons in optical lattices. 
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1. Introduction 

1.1. Bose-Einstein condensation and a new quantum era 

The achievement of Bose-Einstein condensation (BEC) of dilute gases in 1995 [1, 2, 3] 
has marked the beginning of a new era in atomic, molecular and optical physics 
and quantum optics. For the AMO community it was immediately clear that the 
specific experimental techniques of these fields could be used to study problems 
usually encountered in condensed matter physics: degenerate quantum many body 
systems. The condensed matter community remained at this stage much more 
skeptical and arguc^l that at the very end what was achieved experimentally was 
the regime of weakly interacting Bose gases, that had been thoroughly investigated by 
condensed matter theorists in the 50's and 60's [4, 5]. For solid state/condensed matter 
experts the very fact that the AMO experiments dealt with confined systems of finite 
size and typically inhomogeneous density was of technical, rather than fundamental 
importance. Nevertheless, the Nobel foundation decided to give its yearly prize in 
2001 to E. A. Cornell, C. E. Wieman and W. Ketterle "for the achievement of Bose- 
Einstein condensation in dilute gases of alkali atoms, and for early fundamental studies 
of the properties of the condensate" [6, 7]. Today, from the perspective of some years, 
wc sec that due to the efforts of the whole community these fundamental studies 
have enriched amazingly the standard "condensed matter" understanding of static 
and dynamical properties of weakly interacting Bose gases [8] . 

At the same time, the AMO community continued the efforts to extend the BEC 
physics toward new regimes and new challenges. The progress in this directions was 
indeed spectacular and in the beginning of the third millennium is clear both for 
AMO and condensed matter communities that we are entering a truly new quantum 
era with unprecedented possibilities of control on many body systems. In particular 
it became clear that the regime of strongly correlated systems may be reached with 
ultracold atoms and/or molecules. Few years after the first observation of BEC, atomic 
degenerate Fermi gases [9, 10, 11, 12] have been achieved. This has paved the way 
toward the observations of Fermi superfluidity (described in the weak interaction limit 
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by Bardeen-Cooper-Schrieffer (BCS) theory [5]), and the so called BEC-BCS crossover 
in the limit of strong correlations (for recent reviews of an enormous activity in this 
field see [13, 14]). Even earlier, following the seminal proposal by Jaksch et al. [15], 
Greiner et al. [16] observed the signatures of the quantum phase transition from the 
superfluid to the so-called Mott insulator state for bosons confined in an optical lattice. 

Nowadays, ultracold atomic and molecular systems are at the frontier of modern 
quantum physics, and are seriously considered as ones that offer more control than 
solid state systems. It is generally believed that these systems will find highly 
nontrivial applications in quantum information (either as quantum simulators, i.e. 
quantum computers for a special purpose, or as universal ones) or quantum metrology. 
At the level of theory a fascinating "grand unification" takes place: AMO, condensed 
matter, nuclear physics, and even high energy physics theorists join their efforts to 
work on ultracold gases (for recent reviews see [13, 14, 17]). 

1.2. Interactions 

Although quantum gases are very dilute systems (with densities typically ranging 
from 10^^ to 10^^ cm~^), most of their properties are governed by the interaction 
between particles. Usually, in the ultracold regime characteristic of quantum gases 
(temperatures in the nanoKelvin range), only s-wave scattering between particles can 
take place. This allows one to replace the real interatomic potential (which at long 
distances is the usual van dcr Waals interaction) by a pseudo-potential, which is short 
range, isotropic, and characterized by a single parameter, the s-wave scattering length 
a. This contact interaction potential reads 

C/contact(r) = ^^5{r) ^ gSir), (1.1) 

where m is the atomic mass. For a large number of atomic species, the magnitude 
and even the sign of the scattering length a can be tuned by means of an external 
magnetic field. This phenomenon called Feshbach resonance [18, 19] has found multiple 
applications since its first experimental observation in cold gases [20, 21]. 

Despite its simplicity, the interaction potential (1.1) is responsible for an 
extremely rich variety of physical properties of quantum gases. As we mentioned 
above, already for weakly interacting Bose gases, well described by mean-field theory, 
the interactions play a crucial role in the static and dynamic properties of Bose- 
Einstein condensates [8, 22]; one of the most fascinating properties they are responsible 
for is the superfluid character of those gases. The interactions between particles play 
obviously an even more crucial role in the very active field concerning the study of 
strongly correlated systems realized with ultracold atoms [14, 17]. 

For all those reasons, in the last few years, there has been a quest for realizing 
quantum gases with different, richer interactions, in order to obtain even more 
interesting properties. Several researchers, among which| K. Rz§zewski, G. V. 
Shlyapnikov, P. Zoller, G. Kurizki, L. You, D. DeMille, M. A. Baranov, P. Meystre, 
H. Pu, and some of us, have pointed out that the dipole-dipole interaction, acting 
between particles having a permanent electric or magnetic dipole moment, should 
lead to a novel kind of degenerate quantum gases already in the weakly interacting 
limit. Its effects should be even more pronounced in the strongly correlated regime. 

The dipole-dipole interaction has attracted a huge interest for two reasons: 



X We apologize in advance if we have forgotten someone. 
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• Significant experimental progress was made in recent years in the cooling and 

trapping of polar molecules [23], and of atomic species having a large magnetic 
moment. For the case of polar molecules, a very promising technique is to 
associate ultracold atoms by means of Feshbach resonances, and then to use 
photoassociation to bring the weakly bound Feshbach molecules to their ground 
state [24]. In practice this technique requires a very good phase stability of the 
involved lasers. A few months ago D. Jin and J. Ye groups at JILA have been 
able to create a gas of motionally ultracold Rubidium-Potassium molecules in 
their ground rotovibrational state [25]; similar work was done in the group of 
M. Weidemiiller with LiCs molecules [26]. These amazing achievements open 
the way toward degenerate gases with dominant dipole-dipole interactions. For 
the case of magnetic dipoles, Bose-Einstein condensation of ^^Cr, a species with 
a large magnetic moment of 6/iB, was achieved in 2004 [27], and has since then 
allowed for the first experimental investigations of the unique properties of dipolar 
quantum gases [28]. Although the relative effect of the dipole forces in Chromium 
can be tuned using the Feshbach resonance technique, they are typically smaller, 
or at most on the same order as the van der Waals forces. Nevertheless their 
influence on the physics of the Chromium BEC is stunning, as we shall see in the 
following. 

• The properties of the dipole-dipole interaction are radically different from the 
ones of the contact interaction [29]. Indeed, one directly sees from expression 
(2.2) below, giving the interaction energy between two dipoles polarized along the 
same direction, that the dipole-dipole interaction is long-range (it decays like 
where r is the distance between the particles) , and anisotropic (the strength and 
sign of the interaction depend on the angle 9 between the polarization direction 
and the relative position of the particles). Note that if one limits oneself to 
neutral particles, the dipole-dipole interaction is the only interaction between 
electric or magnetic multipole moments which is long-range (interactions between 
higher order multipoles decay fast enough at large distances so that they can be 
replaced by a short range contact pseudo-potential at low temperatures). Long 
range and/or anisotropic interactions are known, already in classical fluids, to lead 
to completely new physical phenomena (sec for example the case of ferrofluids [30] 
in figure 1). Similarly, anisotropy of interactions lies behind the fascinating 
physics of liquid crystals [31]. As we will argue in this review, dipole interactions 
in quantum gases lead also to a variety of novel, fascinating, and sometimes 
completely unexpected effects. 

This review is organized as follows. After a brief description of the main properties 
of the dipole-dipole interaction (section 2), and of the systems in which it can be 
studied (section 3), the mean field theory for a weakly interacting, polarized dipolar 
condensate is presented (section 4). We derive here the non-local Gross-Pitaevski 
equation and discuss its applicability. In subsequent sections we describe a number of 
properties of dipolar BECs, in particular their static and dynamic properties (sections 
5 and 6). In section 7 we enter the very rich field of non-linear, non-local atom optics 
with dipolar gases. Section 8 is devoted to the physics of dipolar spinor condensates. 
Finally, strongly correlated systems obtained by loading a dipolar BEC into an optical 
lattice are described in section 9. Because of the lack of space, some very interesting 
topics are not addressed here; in particular, for a review of the properties of dipolar 
Fermi gases, and of (strongly correlated) rapidly rotating dipolar condensates, the 
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Figure 1. The Rosensweig instability [30] of a fcrrofluid (a colloidal dispersion 
in a carrier liquid of subdomain ferromagnetic particles, with typical dimensions 
of 10 nm) in a magnetic field perpendicular to its surface is a fascinating example 
of the novel physical phenomena appearing in classical physics due to long range, 
anisotropic interactions. Figure reprinted with permission from [32], Copyright 
2007 by the American Physical Society. 



reader is referred to the recent review article [33] . 



2. Dipole-dipole interaction 

2.1. Properties of the dipole-dipole interaction 

For two particles 1 and 2 with dipole moments along the unit vectors ei and 62, 
and whose relative position is r (see figure fa), the energy due to the dipole-dipole 
interaction reads 

rr / N C'dd (ei ■e2)r^ -3(ei -r) (ca ■ r) 

Udd[r) = g . (2.f) 

47r r° 

The coupling constant Cdd is A*oM^ foi' particles having a permanent magnetic dipole 
moment /i (/ig is the permeability of vacuum) and /eq for particles having a 
permanent electric dipole moment d (eq is the permittivity of vacuum). For a polarized 
sample where all dipoles point in the same direction z (figure fb), this expression 
simplifies to 

TT (^\ C'dd 1 - 3 cos^ e 

Udd[r) = , (2.2) 

47r r-^ 

where is the angle between the direction of polarization and the relative position 
of the particles. Two main properties of the dipole-dipole interaction, namely its 
long-range {'-^ ^/t^) and anisotropic character, are obvious from (2.f) and (2.2), and 
contrast strongly with the short-range, isotropic contact interaction (l.f) usually at 
work between particles in ultracold atom clouds. 



Long-range character — In a system of particles interacting via short-range 
interactions, the energy is extensive in the thermodynamic limit. On the contrary, 
in systems with long-range interactions, the energy per particle does not depend only 
on the density, but also on the total number of particles, ft is easy to see that 
a necessary condition for obtaining an extensive energy is that the integral of the 
interaction potential U (r) 

/"OO 

/ f/(r-)d^r, (2.3) 

where D is the dimensionality of the system and rp some short-distance cutoff, 
converges at large distances. For interactions decaying at large distances as l/r", 
this implies that one needs to have i3 < n in order to consider the interaction to be 



CONTENTS 



7 




Figure 2. Two particles interacting via the dipole-dipole interaction, (a) Non- 
polarized case; (b) Polarized case; (c) Two polarized dipoles side by side repel each 
other (black arrows); (d) Two polarized dipoles in a 'head to tail' configuration 
attract each other (black arrows). 

short-range. Therefore, the dipole-dipole interaction (n = 3) is long-range in three 
dimensions, and short range in one and two dimensions. For a more detailed discussion, 
including alternative definitions of the long-range character of a potential, the reader 
is referred to [34]. 

Anisotropy — The dipole-dipole interaction has the angular symmetry of the 
Legendre polynomial of second order P2{cos9), i.e. d-wave. As 9 varies between 
and 7r/2, the factor 1 — 3cos^ 6 varies between —2 and 1, and thus the dipole-dipole 
interaction is repulsive for particles sitting side by side, while it is attractive (with 
twice the strength of the previous case) for dipoles in a 'head-to-tail' configuration (see 
figure 2(c) and (d)). For the special value 6'„i = arccos (1/V3) ~ 54.7° — the so-called 
'magic-angle' used in high resolution solid-state nuclear magnetic resonance [35, 36] — , 
the dipole-dipole interaction vanishes. 

Scattering properties — Usually, the interaction potential between two atoms 
separated by a distance r behaves like —Cg/r^ at large distances. For such a van 
der Waals potential, one can show that in the limit of a vanishing collision energy, 
only the s-wave scattering plays a role. This conies from the general result stating 
that for a central potential falling off at large distances like 1/r", the scattering phase 
shifts 6i{k) scale, for fc ^ 0, like fc^^+i if ^ < (n - 3)/2, and like fc"-^ otherwise [37]. 
In the ultracold regime, the scattering is thus fully characterized by the scattering 
length a. In the study of quantum gases, the true interaction potential between the 
atoms can then be replaced by a pseudo-potential having the same scattering length, 
the so-called contact interaction given by (1.1). 

In the case of the dipole-dipole interaction, the slow decay as 1/r^ at large 
distances implies that for all i, 6e ^ k at low momentum, and all partial waves 
contribute to the scattering amplitude. Moreover, due to the anisotropy of the dipole- 
dipole interaction, partial waves with different angular momenta couple with each 
other. Therefore, one cannot replace the true potential by a short-range, isotropic 
contact interaction. This specificity of the dipolar interaction has an interesting 
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Figure 3. Tuning of the dipole-dipole interaction can be obtained by making the 
dipoles precess around z using a rotating field. 



consequence in the case of a polarized Fermi gas: contrary to the case of a short- 
range interaction, which freezes out at low temperature, the collision cross section for 
identical fermions interacting via the dipole-dipole interaction does not vanish even 
at zero temperature. This could be used to perform evaporative cooling of polarized 
fermions, without the need for sympathetic cooling via a bosonic species. 

Dipolar interactions also play an important role in determining inelastic scattering 
properties. In particular, because of its anisotropy, the dipole-dipole interaction can 
induce spin-flips, leading to dipolar relaxation. The cross-section for dipolar relaxation 
scales with the cube of the dipolc moment [38] , and therefore plays a crucial role in 
strongly dipolar systems (see section 3.4.1). Dipolar relaxation is usually a nuisance, 
but can in fact be used to implement novel cooling schemes inspired by adiabatic 
demagnetization as described in section 3.4.3. 

Fourier transform — In view of studying the elementary excitations in a dipolar 
condensate, as well as for numerical calculations, it is convenient to use the Fourier 
transform of the dipole-dipole interaction. The Fourier transform 



where a is the angle between k and the polarization direction (see Appendix A). 
Remarkably, in three dimensions, the Fourier transform of the dipole-dipole interaction 
does not depend on the modulus of the wavevector k, a feature which is shared by the 
contact interaction (1.1), whose Fourier transform is simply g. 

2.2. Tuning of the dipole-dipole interaction 

By using a rotating polarizing field, it is possible, by time-averaging, to tune the 
dipole-dipole interaction, namely to reduce its effective strength and even change 
its sign [39]. For definitencss we consider here the case of magnetic dipoles /i. in 
a magnetic field B{t) = Be{t) (see figure 3). The unit vector e{t) — costpez + 




(2.4) 



of (2.2) reads 




(2.5) 
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sirup [cos{^t)ex + sm{U,t)ey] is rotated about the 2;-axis on a cone of aperture 2^p at 
an angular frequency 17 which is smaU compared to the Larmor frequency iiB/h, 
but much larger than the trapping frequencies. Then, only the time-average over a 
period 27r/0 of the dipole-dipole interaction (2.1) with ei = 62 = e{t) plays a role to 
determine the properties of the gas. This time-averaged potential reads 



Cdd l-3cos2 6' 



3 cos^ (f — 1 



(2.6) 



47r 

The last factor between brackets decreases from 1 to —1/2 when the tilt angle varies 
from to 7r/2, and vanishes when cp is equal to the magic angle The 'inverted' 
configuration {ip > 6^) in which the averaged dipole-dipole interaction is attractive 
for particles sitting side by side, allows to explore otherwise inaccessible physics (see 
section 7 for some examples of applications). 

3. Creation of a dipolar gas 

In order to realize a quantum gas with significant dipole-dipole interactions, one can 
use particles having either an electric dipole moment rf, or a magnetic dipole moment 
jji. Usually, the dipolar coupling is much higher in the electric case. Indeed, the typical 
order of magnitude of d for an atomic or molecular system is d ~ QeOO) where qf. is the 
electron charge and ao the Bohr radius, while the magnetic moments are on the order 
of the Bohr magneton /ie. Using the definitions of oo and /ie in terms of hmdamental 
constants, one sees that the ratio of magnetic to electric dipolar coupling constants is 

|^~„=~10- (3.) 

where a ~ 1/137 is the fine structure constant. 

For a given species, it is convenient to define various quantities to quantify the 
strength of the dipolar interaction. Prom the dipole moment {i.e. the dipolar coupling 
constant Cdd) and the mass m of the particle, one can define the following length: 

add=Y^. (3.2) 

This 'dipolar length' is a measure of the absolute strength of the dipole-dipole 
interaction. However, in some circumstances, it is the ratio 

add C'dd /o o\ 

£dd = = ^ (3.3) 

of the dipolar length to the s-wave scattering length, comparing the relative strength 
of the dipolar and contact interactions, which determines the physical properties of 
the system. This dipolar parameter needs to be non negligible if one wants to observe 
dipolar effects. The numerical factors in (3.2) are chosen in such a way that for £dd ^ 1 
a homogeneous condensate is unstable against 3D collapse (see section 5.1). Table 1 
summarizes some typical numerical values of the dipolar constants for various atomic 
and molecular species. 

In this section, we review the different systems that can be used in principle to 
study experimentally the effect of the dipole-dipole interaction in degenerate quantum 
gases. We first address the various candidates having an electric dipole moment, 
either static or induced by a laser. The case of magnetic dipoles (the only system 
to date in which strong dipolar effects in a quantum gas have been observed) is then 
described, with an emphasis on the experimental techniques used to achieve Bose- 
Einstein condensation of Chromium. 
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Table 1. Dipolar constants for various atomic and molecular species. For the 
molecular species, the (yet unknown) scattering length is assumed to be 100 ao 
(as the Cs coefficient of the dimer is comparable to the one of a single atom, the 
order of magnitude of the scattering length is similar, but obviously the actual 
value highly depends on the details of the potential) . 



Species 


Dipole moment 






87Rb 




0.7ao 


0.007 


52 Cr 


e.OAtB 


16 ao 


0.16 


KRb 


0.6 D 


2.0 X lO-^o 


20 


NDs 


1.5 D 


3.6 X 10-^ ao 


36 


HON 


3.0 D 


2.4 X 10* ao 


240 



3.1. Polar molecules 

Due to their strong electric dipole moment, polar molecules arc ideal candidates to 
show dipolar effects. Three requirements need to be fulfilled in order for a molecule 
to have a significant dipole moment: 

(i) a heteronuclear molecule, having a permanent dipole moment, is needed§; 

(ii) the molecule must be in a low rovibrational state in order to have a dipole moment 
whose magnitude is not vanishingly small (as would be the case for a highly excited 
vibrational state, especially for Feshbach molecules; indeed, the dipole moment 
scales asymptotically as with the internuclear separation R [42].) and to be 
stable against coUisional relaxation; 

(iii) an external electric field (with a typical value on the order of 10^ V/cm) must 
be applied to orient the molecule in the laboratory frame and approach the 
asymptotic value of the permanent electric dipole moment along the internuclear 
axis (indeed, the ground state J = is rotationally symmetric and therefore 
the dipole moment averages to zero; only via a mixing with higher rotational 
levels, induced by the electric field, does the average dipole become non-zero, see 
Appendix B). Note that this effect can be used to tune the strength of the dipole- 
dipole interaction (but not its sign, unlike the rotating field method described in 
section 2.2). Using additional microwave fields allows for advanced tailoring of 
the interactions between molecules [43]. 

If these requirements are met, the dipole moment is on the order of one Debye 
(1 D ~ 3.335 X 10"^'^ C • m). Assuming that the order of magnitude for the scattering 
length is similar to that of atoms commonly used in BEC experiments (typically 
around 100 oq), the corresponding value of £dd is on the order of 100 (see Table 1), 
meaning that the properties of such a quantum gas would be dominated by the dipole- 
dipole interaction. 

Quantum degenerate gases of polar molecules are a 'Holy Grail' of experimental 
molecular physics. Progress has been made recently in cooling of molecules, with 
techniques such as Stark deceleration (see e.g. [44] for a review) or buffer-gas 
cooling [45, 46, 47], but the densities and temperatures achieved so far are still orders 
of magnitude away from the quantum degenerate regime. A very promising approach 
to degeneracy, actively explored by several groups [24, 25, 26] is to start from already 

§ Note however that exotic homonuclear molecules, such as the ultra-long-range Rydberg molecules 
predicted in [40] and recently observed in [41] can have a permanent dipole moment. 
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ultracold atomic mixtures, and then use a Feshbach resonance to create heteronuclear 
molecules [18]. Created in a highly excited vibrational state, they must then be 
brought to the vibrational ground state, e.g. by photoassociation using STIRAP 
processes, as demonstrated recently [25, 26]. 

3.2. Rydherg atoms 

Extraordinarily large electric dipole moments can be obtained for highly excited 
Rydberg atoms. As the Kepler radius — and thus the dipole moment — scales with n^, 
where n is the main quantum number, the dipolar interaction energy can in principle 
scale like n!^. Individual Rydberg atoms experience lifetimes which scale with n~^. 
However due to the weak binding of the valence electrons and the strong and partially 
attractive forces between Rydberg atoms, the lifetime of a dense gas is limited to time 
scales much shorter than the lifetime of a free Rydberg atom [48] . Therefore Rydberg 
atoms in a BEC [49] are currently investigated as a frozen gas. Collective behaviour in 
the excitation dynamics has been observed, as well as the excitation blockade due to 
dipolar interactions [50] . However hydrodynamic collective phenomena due to moving 
dipoles have not been observed to date. Besides the static dipolar interaction also van 
der Waals interactions (cx n^^) and AC dipolar interactions can occur if neighboring 
energy levels allow for resonant energy transfer via a so-called Forster resonance. 

3.3. Light-induced dipoles 

Atoms in their ground state, which is a parity eigenstate, do not possess an 
electric dipole moment. Their electric polarizability is usually very small, such 
that extreme electric field strengths would be necessary to induce a sizable dipolar 
interaction [51, 52, 53]. Following G. Kurizki and coworkers, one might consider to 
use resonant excitation of a dipole optical allowed transition to induce an AC dipole 
moment on the order of one atomic unit coq. However as this dipole moment also 
couples to the vacuum modes of the radiation field, the spontaneous light forces scale 
just like the light induced dipolar interactions which makes their observation very 
difficult. Nevertheless, the anisotropic nature of the interaction might be used for 
a proof of principle experiment, which would allow to discriminate the spontaneous 
light forces from the dipolar forces [54]. Such interactions have the same form as 
retarded interactions between two dipoles [55]: they contain 1/r^ and radiative 

1/r terms multiplied by the appropriate factors oscillating with the spatial period of 
the laser wavelength. Using an arrangement of several laser fields it has been proposed 
to cancel all anisotropic 1/r^ terms, leaving an effective isotropic, gravity- like 1/r 
potential [56, 57, 58, 59]. In some situations this may lead to self-trapping of the BEC. 
Even before the discovery of the roton instability [60] discussed in section 5.6, a similar 
effect was predicted in a gas with laser induced dipole-dipole interactions [57]. Such 
interactions lead naturally to density modulations of the BEC as in supersolid [59], 
and other effects, such as one dimensional compression of the condensate [58], or 
squeezing [61]. Laser-induced interactions lead in particular to interesting density 
modulations in the condensate, somewhat analogous to self-assembled "supersolid" 
(see section 9.4). Due to the above mentioned limitations caused by spontaneous 
emission, these proposals have not been realized yet. The situation in this respect 
might be more promising if one uses CO2 lasers [62]. 
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Figure 4. Scheme of the energy levels of ^ Cr relevant for the realization of a 
^^Cr BEC. The inset (in gray) gives details of the optical pumping. 

3.4- Magnetic dipoles 

In alkali atoms, the maximum magnetic moment in the ground state is of one Bohr 
magneton {(J-b), said thus the magnetic dipolar effects are very weak. However, very 
recently, dipolar effects have been observed in spinor *^Rb condensates (see section 
8 below) and in condensates with a very small scattering length (obtained using a 
Feshbach resonance), either by studying the size of the condensate (case of ^Li, see 
section 5.2 below) or by using atom interferometry (case of ^^K, see section 9.2 below). 

Some other atoms, like Chromium, Erbium, Europium, Dysprosium, and others, 
have a large magnetic moment of several Bohr magnetons in their ground state, and 
thus experience significant magnetic dipole-dipole interaction. Among them, only ^^Cr 
has been Bose-condensed to date [27, 63]. Chromium has a magnetic dipole moment 
of 6/iB, and a scattering length of about 100 ao [64]. This gives £dd — 0.16 [65], which 
allows to observe a perturbative effect of the dipolar interaction on the expansion 
dynamics of the cloud [66]. Here we describe briefly the main steps leading to the 
creation of a ^^Cr BEC, with a special emphasis on the specificities arising from the 
dipole-dipole interaction. 

3.4-1- Creation of a BEC o/^^^Ct Chromium can be laser cooled using the ^83 ^ 
transition at 425.6 nm (see figure 4). However, strong excited state collisions limit the 
density in a magneto-optical trap (MOT) to relatively low values [67]. Therefore, in a 
typical MOT configuration, the steady state atom number in a Cr MOT is limited to 
a few 10^. In addition the cooling transition is not closed, as there is a decay channel 
from the state to the metastable state (via an intercombination transition) 
with a branching ratio of approximately 1 : 250, 000. These facts seem to rule out any 
hope to achieve Bose condensation of Chromium by standard methods. However, the 
atoms in the metastable state have a magnetic dipole moment of 6 /ib which is strong 
enough so that they remain magnetically trapped in the inhomogeneous magnetic 
field configuration of the MOT. One can thus accumulate large atom numbers in the 
metastable state (where they are decoupled from the MOT light and thus do not 
suffer from excited state collisions), and then, at the end of the MOT phase, repump 
them using light at 663.2 nm. In this way, one ends up with more than 10* ground 



CONTENTS 



13 



state atoms magnetically trapped. In [68], the magnetic field configuration of the 
MOT was modified to the one of a lofle-Pritchard trap, allowing a continuous loading 
of a magnetic trap, in which rf evaporative cooling could be performed. 

However, when the density in the magnetic trap becomes too high, one cannot 
gain anymore in phase-space density due to increasing dipolar relaxation. This two- 
body loss mechanism, in which the spin of one of the colliding atoms is flipped (this 
is allowed as the dipole-dipole interaction does not conserve the spin, but only the 
total angular momentum), is especially important for Chromium, compared to the 
case of alkalis, as its cross section scales as the cube of the magnetic dipole moment. 
Typical relaxation rates of /? ~ 10~^^ cm~^/s were measured in ^^Cr at magnetic 
fields of about 1 G [38] , thus preventing the achievement of BEC in the (magnetically 
trapped) low- field seeking state ms = +3. 

The way to circumvent dipolar relaxation is to optically pump the atoms into 
the absolute ground state ms = —3 (via the ^83 ^ ^Pa transition at 427.6 nm, see 
figure 4) and hold them in an optical dipole trap. Then, in the presence of a magnetic 
field such that the Zeeman splitting between adjacent spin states is much higher 
than the thermal energy, dipolar relaxation is energetically forbidden. One can then 
perform evaporative cooling in the dipole trap and obtain a Chromium condensate [27] . 
Recently, an alternative method has been used to obtain ^^Cr condensates, using direct 
loading of the optical dipole trap [63] . 

3.4-2. Feshbach resonances in ^^Cr A very appealing feature of ^^Cr is the existence 
of several Feshbach resonances. These allow to tune the scattering length a, which, 
close to resonance, varies with the applied external magnetic field B as 



where abg is the background scattering length, Bq is the resonance position (where a 

diverges) and A the resonance width. In order to study the effect of the dipole-dipole 
interaction in a BEC, it is of course interesting to use the Feshbach resonance to reduce 
the scattering length towards zero by approaching Bq + A from above, thus enhancing 

£dd- 

In ^^Cr, for magnetic fields B below 600 G, a total of fourteen resonances were 
found by monitoring inelastic losses in a thermal cloud of atoms in the | ''S3, ms = —3) 
state [69]. An accurate assignment of the resonances was possible by considering 
the selections rules and the shifts of the resonances imposed by the dipole-dipole 
interaction only. In contrast to other atomic species, the dipolar contribution is 
therefore dominant as compared to other coupling mechanisms, like second order 
spin orbit coupling, which have the same symmetry. The inclusion of dipole-dipole 
interaction in multichannel calculations [69] gave a theoretical understanding of the 
width A of the various resonances, which turn out to be relatively small (the broadest 
one, located at Bq = 589 G, having a predicted width of A = 1.7 G only). 

In [70], this resonance was used to enhance dipolar effects in a BEC. An active 
control of the magnetic field at the level of 3 x 10~^ in relative value was implemented, 
allowing for a control of a at the level of ~ oq. Figure 5 shows the measured variation 
of o, inferred from the released energy during expansion (see section 6.2). 

3.4-3. Demagnetization cooling The large magnetic dipole moment of ^^Cr is 
responsible for strong spin-flip collisions, which, as we have seen above, prevent 
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Figure 5. (a) Measured variation of the scattering length across the 589 G 
Feshbach resonance in ^^Cr [70]. (b) Absorption images of the condensate after 
expansion for different values of the magnetic field. One clearly observes a decrease 
in size and an increase in ellipticity when a decreases. 



condensation of Cr in a magnetic trap. However, these inelastic collisions can 
be used to implement a novel cooling scheme, proposed in [71] and demonstrated 
experimentally in [72]. This technique is inspired from the well-known adiabatic 
demagnetization used in solid state physics to cool paramagnetic salts [73]. In the 
context of cold atoms this scheme has been proposed for the first time in [74], and 
termed "elevator cooling" . Particularly important was the analysis of the limitations 
of the scheme due to reabsorption effects in the Raman repumping process. 

The principle of this novel cooling mechanism is represented schematically in 
figure 6. Dipolar relaxation introduces a coupling between spin and external degrees 
of freedom. It can thus be used to cool an atomic cloud by letting a sample, initially 
polarized in the lowest energy state (in a field Bq ^ fceTb/Mj where Tq is the cloud 
temperature), relax towards full thermal equilibrium at a field Bi k^T^/ ^l: energy is 
then absorbed by the spin reservoir at the expense of the kinetic energy (see figure 6). 
The temperature of the sample thus decreases, by an amount which can be up to a 
few tens of percents. By optical pumping, the sample can be polarized again, and a 
new cycle can begin. One can also use a continuous cooling scheme, with the optical 
pumping light always on, and a ramp in magnetic field. This scheme can be seen as 
an evaporation process (where one selects the most energetic particles in the cloud) 
in which the evaporated particles have their energy decreased and are then 'recycled' 
by injecting them back into the trap; it is therefore lossless. 

Note that this scheme is applicable for all dipolar species with a large enough 
dipolar relaxation rate. This rate scales as the third power of the electronic spin. 
There could be even a variant of this cooling technique for electric dipole moments 
in heteronuclear molecules. As under optimized conditions any scattered photon in 
the cooling cycle takes more energy than the mean motional energy the number of 
required photons for a certain cooling rate is much lower than in regular laser cooling 
techniques. Therefore the requirements for the closedness of an optical transition used 
here are much less stringent as compared to regular laser cooling techniques. 




Figure 6. Principle of demagnetization cooling, (a) If the magnetic field is large 
enough so that the Zeeman splitting AE between adjacent Zeeman sublevels is 
much higher than the thermal energy k^T, dipolar relaxation is energetically 
suppressed, and the system, fully polarized in mg = —S, is stable, (b) If one 
reduces the magnetic field so that AE ~ k^T, dipolar relaxation occurs, and 
some kinetic energy is converted into Zeeman energy, (c) By applying an optical 
pumping pulse of a~ polarized light, one can polarize the cloud again, but with a 
decrease of the temperature since the optical pumping process deposits an energy 
which is only on the order of the recoil energy. The excess Zeeman energy is taken 
away by the spontaneously emitted photons. 



This scheme has been successfully applied to Cr, allowing for a reduction of the 
cloud temperature by a factor of two (from 20 to 11 /iK), with almost no atom loss [72]. 
This cooling technique is therefore much more efficient than evaporative cooling, where 
the decrease in atom number is large. An important figure of merit for cooling schemes 
in view of obtaining quantum degeneracy is the gain x in phase-space density p per 
atom loss: 

_ dlnp 

^ = -dh^- (^-^^ 
For evaporative cooling, x is limited in practice|| to values about 4. In [72], the 
measured efficiency of demagnetization cooling reached % ~ 11. 

The practical limitations in view of achieving lower temperatures lie essentially 
in the control of the polarization of the optical pumping light, as any residual cr"*" 
component yields a heating of the cloud, and in the control of stray magnetic fields 
at the milligauss level. However, the recoil temperature should be attainable in 
principle with this technique [74], which could be used, in the future, to realize dipolar 
condensates with large atom numbers. Note that the dipolar coupling mechanism 
between spin and motional degrees of freedom demonstrated in this cooling experiment 
is the same as the one employed in the proposals to observe the quantum version of 
the Einstein-de Haas effect, as explained in section 8. 



4. Non-local Gross-Pitaevskii equation 

Dipolar interactions are expected to change many of the properties of the gas, 
even in the non-degenerate case, where thermodynamical quantities can be affected. 

II Using a higher evaporation threshold increases Xi but the evaporation time then increases 
prohibitively. 
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For example, dipolar interactions lead to a shift of the critical temperature for 
condensation [75, 76], which, although negligible for Cr condensates, could be 
significant for strongly dipolar systems made out of polar molecules. However, the 
most dramatic effects of the dipolar interactions arise for pure condensates. In the 
sections below, unless otherwise stated, we thus consider the case of a gas at zero 
temperature. 

4.1. Pseudo-potential and Gross-Pitaevskii equation 

To describe dilute (and therefore weakly interacting) BECs at zero temperature, the 
mean field approach gives extremely good results [8, 22]. In the case of short range van 
der Waals interactions and low energy scattering further simplification can be made, 
namely the van der Waals interaction potential Kdw(T" — f') may be replaced by the 
pseudo-potential 

47rfi2a , d , 

m d\r — r'\ 

This result has been obtained for the first time in the seminal paper by Huang 

and Yang [77] for the case of a gas of hard spheres. It has however more general 
meaning and holds for arbitrary short range potentials. In the language of many- 
body theory it is the result of the T-matrix, or ladder approximation applied to 
many-body systems [4]. It amounts to the rcsummation of diagrams corresponding to 
multiple two-body scattering. Note that when acting on a non-singular function, the 
pseudopotential is not different from the simple contact (Fermi) potential A'Kfi?a5{r — 
r')/m. This is however not true when we deal with singular, yet square integrable 
functions. In fact, strictly speaking, the contact potential A'Kh?a5{r — r')/m is 
mathematically ill-defined [78]. 

In the mean field theory, where the use of contact interactions is legitimate, 
the order parameter ^{r,t) of the condensate is the solution of the Gross-Pitaevskii 
equation (GPE) [8]: 

*fi^ = -|^A7/;+(Kxt+5|V'l')V'- (4.2) 

The non-linear term proportional to g accounts for the effect of interactions within the 
mean-field approximation, and T4xt denotes the external potential. The normalization 
of ip chosen here is J lipl"^ = N, where N is the total atom number. 

4-2. Validity of non-local Gross-Pitaevskii equation 

In the simple man's approach, to include dipolar effects, one just needs to add an extra 
term to the mean-field potential g\'ip\'^ to account for the effect of the dipole-dipole 
interaction, and one gets 

ih-^ = -^AV; + (Fext + 51^1' + $dd) V'- (4.3) 
where $dd is the dipolar contribution to the mean field interaction 

$dd(r, t) = I |^(r', t)f Uadir - r') dV. (4.4) 

This term is non-local (due to the long-range character of the dipolar interaction) and 
makes it much more complicated to solve the GPE, even numerically, as one now faces 
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an integro-differential equation. In the time-independent case, the left-hand side of 
the above equation has to be replaced by yU'^, with the chemical potential, and the 
GPE becomes: 

-|^AV+(yext+5l^r + M^ = MV'- (4.5) 

It should be stressed that, duo to the long range and anisotropic character of the 
dipole-dipole interactions, it is by no means obvious that one can put the pseudo- 
potential and the real potential into one single equation — obviously this implies that 
the long and short range physics can somehow be treated separately. 

Questions concerning the validity of (4.4) were a subject of intensive studies in the 
recent years. In their pioneering papers, L. You and S. Yi [52, 53] have constructed, 
in the spirit of the ladder approximation, a pseudo-potential for the general case of 
anisotropic potentials. Their results were rigorous, but perturbative, in the sense that 
the Born scattering amplitude from the pseudo-potential was reproducing the exact 
one. The conclusion was that, away from shape resonances, the generalized GPE (4.4) 
is valid, and the effective pseudo-potential has the form (assuming for instance that 
we deal with electric dipoles): 

Ves{r - r ) = —5[r - r )— r - r 

m o\r — r'\ 

^ Aneo |r-r'|3 ' ^^'''^ 

where n = {r — r')/\r — r'\, whereas a{d) depends effectively on the strength of the 
dipole moment d = \d\. 

A. Derevianko succeeded to derive a more rigorous version of the pseudo-potential 
with a velocity dependence [79], which was then used in [80] to calibrate the dipole 
interactions. The author predicted that the effects of dipole interactions should be 
significantly enhanced due to the velocity dependent part of the pseudo-potential. 
Unfortunately, these conclusions were too optimistic, due to some incorrect factors in 
the expressions of [79, 81, 82]. 

Further important contributions to these questions came from J. L. Bohn and 
D. Blume [83, 84]. These authors studied the instability and collapse of the trapped 
dipolar gas and compared the mean field (MF) results with diffusive Monte Carlo 
(DMC) calculations. The DMC results agreed quite accurately with the MF ones, 
provided the variation of the s-wave scattering length with the dipole moment was 
properly taken into account. In fact, this dependence had been already noted 
in [52, 53], and can be traced back to the fact that the rigorous form of dipole- 
dipole interactions contains already a contact 6{r) term [55]. Very careful discussion 
of the differences between the GPE approach, and more exact numerical results of 
diffusive Quantum Monte Carlo methods were presented recently by Astrakharchik et 
al. [85]. These authors point out several difference between DMC and GPE results 
in wide range of parameters, especially reflected in the frequency of the low energy 
excitation "breathing" mode. Very recently, D.-W. Wang [86] has managed to derive 
a general effective many body theory for polar molecules in the strongly interacting 
regime. Wang's approach allows to go beyond the Born approximation approach of Yi 
and You. One of the surprising results is that close to shape resonances, anisotropic 
effects of dipole-dipole interactions are strongly reduced. Phonon dispersion relations 
scale as as in the case of a Coulomb gas. 
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4.3. Variational approach and hydrodynamics 

The time independent GPE can be obtained from the minimization of the energy 
functional: 



J3 



r. (4.7) 



In this context the chemical potential appears simply as a Lagrange multiplier 
arising from the constraint on the normalization of the macroscopic wavefunction ip. 
Minimizing the energy functional (4.7) within a space of trial wavefunctions depending 
on a small number of variational parameters is a convenient way to approach in a 
simple manner a variety of problems; a typical example being the trap geometry 
dependence of the stability of a dipolar BEG (see section 5.3). Variational approaches 
can be extended to the time-dependent case by replacing the energy functional (4.7) 
by an appropriate Lagrange action [87, 88]. 

A useful reformulation of the Gross-Pitaevskii equation is obtained by writing 
ip = ^/ncxp(iS), with n the atomic density and S the phase of the order parameter, 
related to the superfluid velocity field hy v = {h/m)VS. Substituting this expression 
in (4.3) and separating real and imaginary parts, one gets the following set of 
hydrodynamic equations: 

|- + V • {nv) = 0, (4.8) 

which is nothing more than the equation of continuity expressing the conservation of 
mass, and an Euler-like equation: 

m- = - V + ,n + 14. + ^dd - ^ ^ j . (4.9) 

The last term in (4.9), proportional to the Laplacian of ^/n is the quantum pressure 
term arising from inhomogeneities in the density and vanishes in the limit of BEGs 
containing a large number of atoms (Thomas-Fermi limit, see section 5.4). 



5. Ground state properties and excitations 



5.1. Homogeneous gas. Phonon instability 

Because of the partially attractive character of the dipole-dipole interaction, the 
stability of a dipolar BEG is a problem that needs to be addressed. Indeed it is 
well known [8] that a homogeneous condensate with attractive contact interactions 
(a < 0) is unstable, as the Bogoliubov excitations have imaginary frequencies at low 
momentum. 

We consider here a homogeneous dipolar condensate, having an equilibrium 
density no. By considering small density and velocity perturbations with frequency w 
and wavevector k, and linearizing the hydrodynamic equations (4.8) and (4.9) around 
equilibrium, one can show that the excitation spectrum is given by 



m 



g + (3 cos^ a — 1) 



+ 



(5.1) 
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Figure 7. (a) A phonon with k perpendicular to the direction of dipoles 
(a = 7r/2) creates planes of higher density (light gray), in which the dipoles 
are in the plane, corresponding to an instability (see section 5.3 for a discussion 
of the geometry dependence of the stability of a trapped dipolar gas), (b) For k 
parallel to the direction of dipoles (a = 0) the dipoles point out of the planes of 
high density; such a perturbation is thus stable. 



which corresponds to the usual Bogohubov spectrum to = ky^gno/m + h?k'^ / (Am-^) [8] 
with the Fourier transform g of the contact interaction (1.1) complemented by the 
one (2.5) of the dipole-dipole interaction. With the definition (3.3) for Cdd, (5.1) 
implies that a dipolar uniform condensate is unstable for £dd > 1, as phonons (fc —^ 0) 
acquire imaginary frequencies, the most unstable situation being in the case of a 
direction of the wavevector perpendicular to the orientation of the dipoles {a = 7r/2). 
At first sight, this might seem counterintuitive: as dipoles side by side repel each 
other, one could conclude (wrongly) that the most unstable phonons correspond to 
those for which k is parallel to the dipoles. Figure 7 shows how one can understand 
intuitively this behaviour. 



5.2. Trapped gas. Elongation of the cloud 

Like in the case of a BEC with contact interactions, in the presence of an external 
trap (usually harmonic in experiments) new properties arise for a dipolar condensate. 
A prominent effect of the dipole-dipole interaction is to elongate the condensate along 
the direction z along which the dipoles are oriented [89, 90, 53]. This magnetostriction 
effect (a change of the shape and volume of the atomic cloud due to internal magnetic 
forces) can be understood in a very simple way for a spherically symmetric trap (of 
angular frequency uj) in the perturbative regime Edd <C 1. To zeroth order, the density 
distribution is given, in the Thomas-Fermi limit, by n{r) — no{l — r^/i?^), where R 
is the Thomas- Fermi radius of the condensate (see figure 8a) . One can then calculate 
to first order in Sdd the mean-field dipolar potential (4.4) created by this distribution; 
one finds [39] 

2 f r2 if r < i? 

$dd(r) -£dd^(l-3cos2 0) <^ i?5 ^ (5.2) 

5 [ if ^ > ^ 

i.e. the dipolar mean field potential has the shape of a saddle, with minima located 
on the z axis (see figure 8b). It is therefore energetically favorable for the cloud to 
become elongated along z. One can actually show that this conclusion remains valid 
even if the cloud is anisotropic, and for larger values of edd [91, 92, 93]. 

Very recently, the spatial extent of a ^Li BEC was studied as a function of the 
scattering length close to a Feshbach resonance [94]. For very small scattering lengths. 
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Figure 8. (a) Inverted parabola density distribution n(r) in tlie Tliomas Fermi 
regime in the absence of dipole-dipole interaction, (b) Saddle-like mean field 
dipolar potential (5.2) induced by the density distribution displayed in (a). 



the elongation effect due to the dipole-dipole interaction could be seen unambiguously, 
in spite of the small value of the magnetic dipole moment. 

5.3. Trapped gas. Geometrical stabilization 

A BEC with pure contact attractive interactions (a < 0) is unstable in the 
homogeneous case, but, in a trap, stabilization by the quantum pressure can occur for 
small atom numbers, namely if 

ffl < 0.58, (5.3) 

dho 



where N is the atom number and Oho = y^h/{muj) is the harmonic oscillator length 
corresponding to the trap frequency oj [95]. Here the trap has been supposed isotropic, 
but, for anisotropic traps, the dependence on the trap geometry is weak [96]. 

The situation is radically different in the case of a BEC with dipolar interactions. 
Due to the anisotropy of the dipole-dipole interaction, the partially attractive character 
of the interaction can be "hidden" by confining the atoms more strongly in the 
direction along which the dipoles are aligned. Let us consider for simplicity a 
cylindrically symmetric trap, with a symmetry axis z coinciding with the orientation 
of the dipoles. The axial (resp. radial) trapping frequency is denoted lu^ (resp. ujp). 
It is then intuitively clear that for a prolate trap (aspect ratio A = uj^/ujp < 1), the 
dipole-dipole interaction is essentially attractive, and in such a trap a dipolar BEC 
should be unstable, even in the presence of a (weak) repulsive contact interaction [see 
figure 9(a)]. On the contrary, in a very oblate trap, the dipole-dipole interaction is 
essentially repulsive, leading to a stable BEC even in the presence of weak attractive 
contact interactions [see figure 9(b)]. One therefore expects that, for a given value 
of A, there exists a critical value flcrit of the scattering length below which a dipolar 
BEC is unstable; from the discussion above, Ocrit should intuitively be a decreasing 
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(a) 



(b) 



Figure 9. Intuitive picture for the geometry-dependent stability of trapped 
a dipolar BEC. (a) In a prolate (cigar-shaped) trap with the dipoles oriented 
along the weak confinement axis of the trap, the main effect of the dipole-dipole 
interaction is attractive, which leads to an instability of the condensate, (b) 
In an oblate (pancake-shaped) trap with the dipoles oriented along the strong 
confinement axis, the dipole-dipole interaction is essentially repulsive, and the 
BEC is stable. 



function of A, and the asymptotic value of Ocrit for A — )■ (resp. A oo) should be 
positive (resp. negative). 

A simple way to go beyond this qualitative picture and obtain an estimate for 
'2crit(A) is to use a variational method. For this purpose, we assume that the condensate 
wavefunction ip is gaussian, with an axial size cr^ and a radial size Up that we take as 
variational parameters: 



1 



^"ho y'^p "2 



(5.4) 



Here, aho = (mw) is the harmonic oscillator length corresponding to the average 
trap frequency a) = {utpLVzY^^. Inserting Ansatz (5.4) into the energy functional (4.7) 
leads to the following expression for the energy: 

E{ap, (Tz) = -Ekin + i?trap + -Eint, (5.5) 

with the kinetic energy 

the potential energy due to the trap 

£;trap = ^(2.^+AV.^), (5.7) 

and the interaction (contact and dipolar) energy 



Eint = — == — /(k) • 

V27raho <^p'^z \add ) 



(5.8) 



The dipolar contribution in the last part is most easily calculated in momentum space 
as 

Edd^lJ "(^)'^(^')t^dd(r - r') dV dV 

lJdd{k)n^{k)d^k, (5.9) 



2(27r)3 

where n(k) is the Fourier transform of the density distribution (and therefore, in this 
case, still a Gaussian). In (5.8), k — <7p/<7z is the aspect ratio of the cloud (which 
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Figure 10. The function f{x) entering the calculation of the dipolar mean-field 
energy. 



differs from the one of the trap due to the elongation induced by the dipole-dipole 
interaction as discussed above), and / is given by 

2^2 SK^artanh-y/l — 



/(^) 



1 



(5.10) 



1 - k2 (1 - ^2)3/2 

The function /(k), displayed in figure 10, is monotonically decreasing, has asymptotic 
values /(O) = 1 and /(oo) = — 2, and vanishes for k = 1 (implying that for an isotropic 
density distribution the dipole-dipole mean- field potential averages to zero). 

To determine the stability threshold acrit(A), one needs to minimize (5.5) with 
respect to Up and az for fixed values of N, X and a). For a > flcrit, one has a (at 
least local) minimum of energy for finite values of cfp^z, while as soon as a < fldit, no 
such minimum exists. Figure 11 shows contour plots of E{ap,az) for N — 20,000, 
A = 10 and different values of a, clearly showing that acrit(lO) ~ —8.5 ao for the chosen 
parameters. In figure 13, the critical scattering length acrit(A) obtained in this way is 
shown as a thick line for Cu = 2ttx 800 Hz and N — 20, 000 atoms. In the limit N oo, 
the asymptotic behaviour of this curve (0^4^(0) = Odd and 0^1^(00) = — 2add) can be 
easily understood, as only the sign of the interaction term (5.8) (which scales as TV^ 
and not as N like the kinetic and potential energy) determines the stability. For an 
extremely pancake-shaped trap A — > 00, the cloud has an aspect ratio k — s- 00, and, 
as lim f{x) = —2, the condensate is (meta-)stable only if a > — 2add- In the same 

way, one readily understands that for A ^ 0, the critical scattering length is add- The 
minimal value of A for which a purely dipolar condensate (a — 0) is stable is the one 
for which k = 1 and is found numerically to be close to A ~ 5.2 [89, 97, 53, 92, 98]. 

In [98], the influence of the trapping geometry on the stability of a ^^Cr BEC was 
investigated experimentally. A combination of an optical dipole trap and of one site of 
a long period (7 /im) optical lattice provided an harmonic trap cylindrically symmetric 
along the z direction (along which the dipoles are aligned), with an aspect ratio A 
that could be varied over two orders of magnitude (from A ~ 0.1 — prolate trap — 
to A ~ 10 — oblate trap — ), while keeping the average trap frequency a) = (w^oj^)^/'^ 
almost constant (with a value of 27r x 800 Hz). Using the Feshbach resonance at 589 G, 
the scattering length was ramped adiabatically to a final value a and the atom number 
in the BEC was measured. A typical measurement is shown in figure 12. When a 
is reduced, the atom number decreases, first slowly, and then very abruptly when a 
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Figure 11. The energy landscape E{ap,az) as a function of tlie variational 
parameters CTp and of the Gaussian Ansatz, for a trap of aspect ratio A = 10, 
and various values of the scattering length a. When a decreases, one goes from a 
global minimum (at a = 18ao) to a local minimum corresponding to a metastable 
condensate (at a = lOao). This local minimum vanishes at a = Ocrit (here 
— 8.5ao). Below acrit, the energy can be lowered without bound by forming an 
infinitely thin cigar-shaped cloud. 
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Figure 12. Experimental observation of the geometry-dependent stability of a 
dipolar BEG. (a) BEG atom number N as a. function of a for a spherical trap; N 
vanishes for a smaller than a^rit — ISoq. (b) For an oblate trap (A = 10), one 
has Ocrit — — 2ao; such a trap can thus stabilize a purely dipolar BEG. In (a) and 
(b) the solid lines are fits to the empirical threshold law (a — cicrit)''- (c) Sample 
images of the atomic cloud as a function of a for A = 10. 



approaches a critical value Ocrit, below which no condensate can be observed. Figure 
13 shows the measured value of Ocrit as a function of A. One clearly observes that 
for prolate traps, Ocrit is close to Odd, as expected from the discussion above, while 
for the most pancake-shaped trap A = 10 the critical scattering length is close to 
zero: for such a geometry, a purely dipolar condensate is stable. The solid line is the 
stability threshold acrit(A) obtained by the gaussian Ansatz for a number of atoms 
N = 2 X W^, which shows a good agreement with the measurements. Note that for 
the parameters used in the experiment, the critical scattering length for pure contact 
interaction, given by (5.3) would be — 0.3ao for A = 1, which clearly shows that the 
instability is driven here by the dipole-dipole interaction. 




Figure 13. Stability diagram of a dipolar condensate in the plane (A, a). The dots 
with error bars correspond to the experimental data [98]; the thick solid line to 
the threshold acrit(-^) obtained using the Gaussian ansatz (5.4) with N = 20, 000; 
the thin solid line to the numerical solution of the GPE (4.3) [99]. 



To calculate the exact stability threshold, one needs to resort to a numerical 
solution of the GPE (4.3); the result of such a calculation [99] is displayed as a thin 
line on figure 13 and shows a very good agreement with the data. The numerical 
solution reveals, for some values of the parameters (A, a) close to the instability region, 
the appearance of 'biconcave' condensates, where the density has a local minimum in 
the center of the trap [100] . ^ 

5.4- Trapped gas. Thom,as- Fermi regime 

As we shall see in the next sections, a very important approximation in the case of 
dipolar gases is the so-called Thomas-Fermi (TF) limit, in which quantum pressure 
effects are neglected. Amazingly, the TF solutions for the ground state of the trapped 
BEG have the same inverted parabola shape as in the case of contact interactions. 
This has been pointed out for the first time in [60] , where, however, the trapping was 
restricted to the z-direction, while in the other directions the systems was assumed 
to be homogenous. The Thomas- Fermi approach was also used to study fermionic 
dipolar gases (see [101] and references therein). 

The exact solutions of the dipolar BEG hydrodynamics in 3D were presented 
in a series of beautiful papers by O'Dell, Giovanazzi and Eberlein [91, 92]. These 
authors have used oblate spheroidal coordinates and solved the TF equations for the 
ground state in cylindrically symmetric traps. They have also considered the stability 
of the three most relevant perturbations: local density perturbations, "scaling" 

^ The experimental observation of such biconcave condensate (which show in a striking manner the 
long-range character of the dipole-dipole interaction) is difficult for several reasons: (i) the density 
dip does not survive in time of flight (which implies that in-situ imaging would be needed to detect 
it) (ii) it has a small contrast of only a few percent and (iii) the regions in the plane (A, a) where 
the biconcave condensate exist have a very small area. However, the use of potentials flatter than 
harmonic traps, such as a quartic or a box-like potential, should relax considerably the constraints 
(ii) and (iii) (S. Ronen, 2008, private communication). 
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Figure 14. Aspect ratio k of tlie condensate as a function of the dipole-dipole to 
s-wave coupling ratio Sdj. Each line is for a different trap aspect ratio 7 = w^/o^x, 
which can be read off by noting that /t(£dd = 0) = 7- When < k < 1 the 
condensate is prolate; for k > 1 it is oblate. Likewise, for < 7 < 1 the trap is 
prolate, and when 7 > 1 the trap is oblate (figure courtesy of C. Eberlein). 



perturbations, and "Saturn-ring" perturbations. 

In particular, the ground state density in the cyhndrically symmetric case has the 
form: 



n(r)=no(l-|^-J, 



(5.11) 



for n{r) ^ 0, where no = 15N/{8TrR'^Rz). These expressions are exactly the same 
as in the case of contact interactions. The difference is, of course, in the explicit 
expressions for the radii: 

1/5 



Rt. — R 



V 



4:Trmu)?. 



1 + Edd 7; - 



2 l-«;2 



1 



(5.12) 



and Rz = Rx/k. The condensate aspect ratio k is determined by the transcendental 
equation 

^ /(«) 



3K£dd 



+ 1 



- 1 



+ (edd-l)(K'-a;,7a;^) = 0, (5.13) 



where /(k) is defined in (5.10). A plot of the condensate aspect ratio as a function 

of Edd is shown in figure 14. The TF approach is also extremely useful to study the 
dynamics of dipolar condensates, for instance their free expansion (see section 6.1). 



5.5. Trapped gas. Excitations 

The unusual properties of the ground state of dipolar condensates have their 
counterpart in the excitations of the system. They are expected to exhibit novel 
character and symmetries, as well as new types of instabilities. Indeed, as we shall 
see in the next subsection, even in the pancake traps with dipole moments polarized 
orthogonally to the pancake plane, the excitation spectrum reveals an instability, at 
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the so-called roton-maxon minimum. Before discussing the pancake case, let us first 
consider in this section the case of a moderate aspect ratio of the trap. 

The study of the excitations of a dipolar BEC should in principle be realized using 
the non-local Bogoliubov-de Gennes (BdG) equations. Such an approach is however 
technically very difficult, and for this reason approximate methods arc useful. We 
discuss here the results of pioneering papers (see [52, 53, 97] and references therein). 
Goral and Santos [97] apply the dynamical variational principle, developed for the 
contact GPE in [87, 88]. The idea consists in writing the time dependent condensate 
wave function in the Gaussian form 

xP{x,y,z,t)=A{t) exp[-^7V2^«^(^)'-*^'/3r,(^)] , (5.14) 

with time-dependent variational parameters describing the Gaussian widths Wri{t), 
and phases Pn{t), while A(t) takes care for normalization. The dynamical variational 
principle implies equations of motions for the widths and phases. Stationary solutions 
of these equations describe the ground state BEC, small deviations from the ground 
state describe the lowest energy excitations. Goral and Santos consider the case of a 
polarized dipolar gas and discuss the influence of dipole-dipole forces on the stability 
of the condensate and the excitation spectrum. The authors extend their discussion 
of the ground state and excitations properties to the case of a gas composed of the 
two anti-parallel dipolar components. 

One of the most interesting results of this paper is the study of the nature of 
the collapse instability. In the standard case of contact interactions three modes are 
relevant at low energies (see figure 15): two "quadrupole"-like modes (1 and 3), and 
one "monopole" mode (2). It is the latter one which becomes unstable at the collapse 
(when the scattering length changes sign). The frequency of the breathing monopole 
mode 2 goes to zero in a certain manner, namely, if 7 denotes the ratio of the non-linear 
energy to the trap energy, and 7c is correspondingly its critical value, the frequency of 
the breathing monopole mode 2 goes to zero as [7— 7c|^^^ [102, 103] when 7 approaches 
criticality from below. 

In the case of dipolar gases with dominant dipole interaction the situation is 
similar only for aspect ratios A -C 1 far from criticality, where the lowest frequency 
mode is the breathing mode, and its frequency tends to zero as |7dd — 7c,dd|'', 
with /? ^ 1/4, when the energy of the dipolar interactions approach the criticality. 
Numerical analysis indicates that as one approaches the criticality the exponent (3 
remains close to 1/4, but the geometry of the zero frequency mode is completely 
different: it attains the quadrupole character and becomes a superposition of modes 
1 and 3. Very close to the criticality the exponent /3 grows up to the value ~ 2. 
These results imply already that one should expect for the dipolar gas a completely 
different character of the collapse dynamics than for a gas with contact interactions 
and negative scattering length. We will discuss it in detail in the next section. 

The variational method employing the Gaussian ansatz is evidently the simplest 
approach to study the excitations and dynamics of the dipolar BECs, and for this 
reason it was used by many authors. After the seminal experiments of the JILA 
group with *'''Rb, in which efficient, i.e. practically loss-free control of the scattering 
length was achieved [104], Yi and You used this method to investigate the possibility 
of observing dipolar effects in the shape oscillations and expansion of a dipolar 
condensate [105, 106]. 

Direct solution of the Bogoliubov-de Gennes equations is, as we said, difficult, 
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mode 1 mode 2 mode 3 



Figure 15. Graphical representation of the oscillation modes of the condensate; 
the modes 1 and 3 are "quadrupole"-like, the mode 2 is the breathing "monopole" 
mode. 

but not impossible. First of all, they become easier to solve in the Thomas-Fermi 
(TF) limit. The first step toward it was made in the "roton-maxon" paper of Santos 
et al. [60]. These authors considered an infinite pancake trap {i.e. a slab) with dipoles 
oriented along the 2;-direction and in the TF limit. It turns out that the density profile 
in the TF limit is given by an inverted parabola, just like in the case of the standard 
BEC [8]. The BdG equations for excitations close to the critical parameter region, 
where the "rotonization" of the spectrum appears, can be solved analytically in terms 
of series of Gegenbauer polynomials, as we discuss in the next subsection. 

Very recently, in a beautiful paper, Ronen et al. [107] developed an efficient 
method for solving BdG equations for the dipolar BEC with cylindrical symmetry. 
The algorithm is very fast and accurate, and is based on an efficient use of the 
Hankel transform. The authors study excitations in different geometries (from cigar 
to pancake; for typical results of BdG spectra, see figure 16) and in particular they 
calculate for the first time the dipolar condensate depletion for various regimes of 
parameters. 

5.6. Trapped gas. Roton-maxon spectrum 

As we have mentioned, a dipolar gas exhibits two kinds of instabilities. In cigar- 
like traps, when the dipole is oriented along the trap axis, the dipolar interactions 
have an attractive character and the gas collapses. The collapse is similar to the 
case of the contact interactions with negative; scattering length, but has a different 
geometrical nature, and different critical scaling behavior. There is, however, another 
instability mechanism that occurs even in quasi-2D pancake traps for dipoles polarized 
perpendicularly to the trap plane (along the z-axis). In this case, when the dipolar 
interactions are sufficiently strong, the gas, despite the quasi-2D trap geometry, feels 
the 3D nature of the dipolar interactions, i.e. their partially attractive character. 
This so-called "roton-maxon" instability has been discovered in [60], and discussed by 
many authors since then. 

In the original paper [60] the authors considered an infinite pancake trap (slab 
geometry) with dipoles oriented along the z-dircction perpendicular to the trap plane. 
The roton-maxon physics occurs in the TF limit in the 2-direction. The condensate 
was hence assumed to have a TF shape along the 2-axis and constant amplitude 
with respect to x and y coordinates. The 3D Bogoliubov-de Gennes equations 
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Figure 16. Excitation frequencies as a function of the dipolar parameter 
D, = {N— l)mCdd/(47r?l-^aijo,p) for a dipolar BEC in the JILA pancake trap [108] 
(LO^jiOp = \/8), with a zero scattering length. Plotted are modes with m = — 4. 
The three lines that extend to higher Dt, are the variational results of [97] (figure 
courtesy of S. Ronen). 

were then solved. Here, we follow a somewhat simplified effective 2D approach 
of [f09], which nevertheless captures the main physics. The bosons are also polarized 
along the z direction, so that the dipolar interaction in momentum space reads 
Vdd = Cddikl/k"^ — 1/3). We introduce also the ratio /3 = 3edd/(47r). As we shall 
see below the cases of repulsive (attractive) contact interactions with g > {g < 0) 
lead to qualitatively different physics. We first consider the case of positive g. 

In the standard quasi-2D approach, one assumes that the condensate has a 
Gaussian shape along the z direction, with the width determined by the harmonic 
potential. The problem is then projected onto 2D, by integrating over the condensate 
profile. Such an approach is valid when the chemical potential fi <C hu. However, the 
roton-maxon instability occurs outside this regime and hence the standard quasi-2D 
approach caimot be employed to understand it. 

However we can still use an effective quasi-2D approach, where we assume that 
the bosonic wave-function in the z direction is given by the TF profile. In this way 
we may obtain an effectively 2D model with the dipole interactions "averaged" along 
the z direction. This approach corresponds approximately to the approach of [60], 
with the only difference that for the lowest energy branch of excitations we neglect 
their "kinetic" energy, i.e. terms in the Bogoliubov-de Gennes equation involving 
derivatives with respect to z. The quantitative differences between the exact results 
of [60], and the present effective quasi-2D approach amount typically to 10-20 percent 
in the entire regime of /3 < 2 (which is the relevant regime for rotonization [60]). 
Qualitatively, both approaches describe the same physics, and the same mechanism of 
appearance of the instability, namely the momentum dependence of the dipole-dipole 
interactions. 

For the purpose of this review we apply yet another approximation, and use for 
simplicity a Gaussian profile in the z direction with a variationally determined width 
£, which turns out to be of order of the TF radius £ ~ £tf / VS. The dipolar interaction 
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Figure 17. The excitation spectrum Q(k^) as a function of k±. Tiie blue, black 
and red curves correspond to 33^ = 2.0, 3.0, 3.44, respectively (figure courtesy of 
O. Dutta). 



in 2D, after integrating out the z direction, takes the form 



where 



V(fci) = l-^J|fciferfc 



exp [kle'/2] , 



(5.15) 



and k'j_ 



— k±£. The Hamihonian that generated the Bogohubov-de 
Gennes equations reduces then in the absence of the contact interactions, i.e. for 
g = 0, to H = X^fcj^ ^{f^±)bl.^bk^, where 6^^^ and b^^ are Bogohubov quasi-particle 
operators. The excitation spectrum is given, in units of the trap frequency, by 



n^~k^) = ^ + g^^v{k^)~e 



(5.16) 



where the dimensionless interaction strength is defined as (73d = mCddnl/ {A'Kh^). 
The interaction in (5.15) is repulsive for small momenta and attractive in the high 
momentum limit (with a zero-crossing at ki_i ~ 1). Due to this fact the properties 
of the excitation spectrum in (5.16) are very different from that of bosons with 
contact interactions. For any bosonic density, r2(A;j_) exhibits two regimes: (i) phonon 
spectrum for small momenta, and (ii) free particle spectrum for higher momenta. 
For (73(j greater than a certain critical value, f2(fc^) has a minimum (see figure 17) 
at the intermediate momentum regime. Following Landau, the excitations around 
the minimum are called "rotons" and il(fco) is known as the "roton gap" [60]. With 
increasing 53d the "roton" gap decreases and eventually vanishes for a critical particle 
density. When the critical density is exceeded, r2(fco) becomes imaginary and the 
condensate becomes unstable. 

Once more we stress that in the presence of repulsive contact interactions, the 
roton instability occurs in pancake traps in the regime in which the standard quasi- 
2D approximation does not hold. The condensate attains the TF profile in the z 
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direction and the systems starts to experience the 3D (partially attractive) nature of 
the dipolar forces [110]. The situation is however different if the contact interactions 
are attractive {g < 0). In that case the roton- instability appears already in the 
standard quasi-2D regime, when the chemical potential /it <C tuv, when the condensate 
profile is Gaussian with the bare harmonic oscillator width. While in the previous case 
it was the attractive part of the dipolar interactions that led to the roton-instability, in 
this quasi-2D situation with g <0 the situation is different. The dipolar interactions 
(which are in average repulsive in this case) stabilize for fi = {g + 8TTgdd/3)n > 
the phonon instability (which in absence of dipole-dipole interaction leads to the well- 
known collapse for gases with a < 0) [111]. Note that /x > but may be kept 
well below hujz- However due to the momentum dependence of the dipolc-dipole 
interaction (which acquires the form (5.15)), one encounters the roton instability for 
— 3/87r < /3 < /3cr (where /3ci depends on the value of g), whereas for /? > 0cr the 
system is stable (as long as it remains 2D [112]). Note that in this case, a larger 
dipole strength stabilizes (sic!) the gas. The reason is that the roton in the quasi-2D 
scenario is not actually induced by the attractive nature of the dipolar interactions at 
large momenta, but by the fact that g < 0. As a consequence, quasi-2D roton may 
occur for much lower momenta than 1/lz- 

The presence of a roton minimum in the spectrum of elementary excitations 
may be revealed in various ways. On one hand, and using Landau superfluidity 
criterion [113], it is clear that the superfluid critical velocity is reduced in the presence 
of a roton minimum [60]. Another alternative experimental signature of the roton 
could be provided at finite temperatures, where the thermal activation of rotons may 
induce a "halo" effect in time-of- flight images [114]. Finally, we would also like to 
point that (as discussed in section 7.3), the presence of even a shallow roton minimum 
may alter dramatically the response of the system against a periodic driving [115]. 

The question whether "there is a life after the roton instability" was studied by 
many researchers, in the hope to find novel types of stable Bose superfluids, that would 
have supersolid character, i.e. self- assembled density modulations. This is suggested 
by the fact that the instability occurs at a specific value of the momentum, indicating 
instability toward a non-uniform ground state [116, 117]. The ultimate answer to 
the question is negative: the condensate undergoes a sequence of local collapses, as 
shown for the first time in [118, 111, 119]. In [111] it is shown numerically that in 
the mean field theory supersolid states of dipolar BEG are unstable. The authors 
of [120] use a variational ansatz with density wave modulations along the z-dircction 
in a cylindrically symmetric trap, and show that it is not stable for a dipolar gas, due 
to the roton instability. It can, however be stabilized, by allowing for admixture of a 
single component polarized Fermi gas. This idea was followed further in [121], where 
the stability of dipolar bosons-fermions mixture in pancake cylindrically symmetric 
traps was investigated at T = using a variational approach. Fermions- induced 
interactions stabilize the system in such traps and allow for quantum phase transition 
from the Gaussian shape BEG to a supersolid state, characterized by a hexagonal 
density wave. 

Interestingly, while fermions stabilize dipolar bosons leading to novel type of 
states, the opposite is also true: boson mediated interactions between polarized 
fermions may lead to appearance of the "exotic" Fermi superfluids with p-wave, /- 
wave, or even h-wavc pairing [109]. All of these states exhibit topological order, 
and admit non-Abelian anionic excitations [122], that can be used for topologically 
protected quantum information processing [123]. 
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The discovery of the roton instabihty in dipolar gases stimulated the search for 
supcrsolid structures and density modulations. Roncn et al. studied angular and 
radial roton instabilities in purely dipolar BECs in oblate traps [100], and have shown 
that in some situations the condensate attains a biconcave density profile (see the end 
of section 5.3). This results have been generalized to the case of finite temperature 
using the Hartree-Fock-Bogoliubov approach [124]. Such structures are very sensitive 
to changes in the trapping potential; their relation to the roton instabilities for the 
case of a trapped dipolar EEC without and with a vortex has been studied in [125]. 
Dutta and Meystre [118] predicted similar effects in anisotropic traps. 

Very recently, the stability, excitations and roton instabilities have been discussed 
for the case of a dipolar BEC with a vortex [126]. Recent observation of the dipolar 
effects in Bloch oscillations with ^^K [127, 128] with s-wave scattering tuned to zero 
stimulated studies of collective excitations and roton instabilities in multi-layer stacks 
of dipolar condensates. In [114, 112], an enhancement of the roton instability was 
predicted. Note also that, interestingly, a gas of light-induced dipoles (see section 3.3) 
were predicted by Kurizki and coworkers to display roton instabilities [57] . 



6. Dynamics of a dipolar gas 

6.1. Self-similar expansion in the Thomas-Fermi regime 

In most experiments on BECs, all the information is obtained from an absorption 
image of the cloud taken after a period of free expansion ("time of flight"), which 
acts as a "magnifier" allowing to resolve optically the BEC (in a trap, the typical size 
of the BEC is on the order of a few microns, making it difficult to image it in situ 
with a good resolution). It is therefore of great practical importance to describe the 
expansion of a condensate released from a trap. 

In the case of a BEC with contact interactions in a harmonic trap of frequencies uji 
{i = X, y, z) in the Thomas- Fermi limit (see section 5.4 above), a remarkable property 
allows for a very simple description of the free expansion [129, 130]: the in-trap density 
profile of the condensate (which is, in the Thomas-Fermi limit, an inverted parabola) 
is merely reseated upon time of flight. The Thomas- Fermi radii Ri{t) at time t are 
given by 

Ri{t) = Ri{Q)hi{t), (6.1) 

where the scaling parameters bi are solution of the following set of coupled diflFerential 
equations: 

where Wi(0) stands for the trap frequency along the direction i before the trap is 
turned off. 

The underlying reason for the existence of such a scaling solution is the fact that, 
for a parabolic density distribution n, the mean-field term gn is also quadratic in the 
coordinates, yielding only quadratic terms in the hydrodynamic equations (4.8) and 
(4.9). One can show that this property remains valid in the case of dipolar condensate 
in the Thomas-Fermi limit, as the following non-trivial property holds: if the density 
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distribution n is parabolic n = no max {0,1 — ^^xf/R"^), then the mean field potential 
(4.4) due to the dipolc-dipole interaction is a quadratic form in the coordinates^. 

This property can be understood in the following way [92]. Starting from the 
identity 



one can prove that 



$dd(r) = -Cdd ( ^ct>{r) + in(r) ) , (6.4) 



where 

' dV. (6.5) 



47r|r - 

The last equality shows that the "potential" fulfills Poisson's equation A<p = —n. 
From this, one deduces that the most general form of 0, when one has a parabolic 
density distribution n, is a polynomial of order two in the variables (a;^, y^, z^), and 
thus, from (6.4), one deduces that $dd is also quadratic in {x,y,z). The actual 
analytical calculation of the coefficients of this quadratic form, carried out for the 
cylindrically symmetric case in [91, 92] and for the general anisotropic case in [93], is 
far from trivial, especially in the latter case, which involves a generalization of (5.10) 
to two cloud aspect ratios. 

One can then generalize Thomas-Fermi scaling equations (6.2) to the case of 
a condensate with both contact and dipolar interactions. The corresponding set of 
differential equations now reads 

= ui^^ul [^ + £ddFib,,by,h)] ii = x,y,z). (6.6) 

where the function F includes the effect of the dipole-dipole interaction [93] . Solving 
these equations, one can for example study the time evolution, during free expansion, 
of the cloud aspect ratio, to reveal the effects of the dipolar interaction (see section 
below). 



6.2. A quantum ferro fluid 

The superfluid hydrodynamic equations describing e.g. the expansion of a condensate 

are modified by the long-range, anisotropic dipole-dipole interaction. In classical 
fluids, such magnetic interactions modifying the hydrodynamic properties can be 
observed in ferrofluids, which are colloidal suspensions of nanometric ferromagnetic 
particles [30]. In that sense, a dipolar condensate can be called a quantum ferro fluid. 
As seen in the preceding section, a clear and simple way to demonstrate the effect of 
the dipole-dipole interaction is to study the expansion of the condensate when it is 
released from the trap. 

By measuring the aspect ratio of the condensate as a function of the expansion 
time for two different orientations of the dipoles with respect to the trap axes, the 
elongation of the cloud along the magnetization direction could be clearly observed 
in [66, 93]. However, since £dd — 0.16 for '"'^Cr away from Feshbach resonances, the 

+ As shown, in the much more restrictive case of a spherical trap and for £;j;j <g 1, in section 5.2 (see 
(5.2) and figure 8). 
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Figure 18. Free expansion of a dipolar condensate for two different orientations 
of the dipoles with respect to the trap axes. The squares are the experimental 
results; the solid lines are the prediction of the scaling equations (6.6) without 
any adjustable parameter, (a) Perturbative regime e^d = 0.16, the dipole-dipole 
interaction only gives a small departure from the contact interaction prediction 
(data taken from [93]). (b) Perturbative regime (for a different trap geometry) 
£dd — 0.16 far above the Feshbach resonance, the effect of the dipole-dipole 
interaction is similar to the case (a), (c) Closer to the resonance, with Edd = 0.75; 
in that case, the dipole-dipole interaction is strong enough to inhibit the usual 
inversion of ellipticity in time of flight. Data in (b) and (c) arc taken from [70] . 

effect, shown in figure 18(a), is small. By using the 589 G Feshbach resonance in 
order to decrease the scattering length a and thus enhance Edd, a much larger effect 
of the dipole-dipole interaction on the expansion dynamics was observed in [70]. For 
example, by increasing £dd to values close to one, the usual inversion of ellipticity of 
the condensate during the time of flight is inhibited by the the strong dipolar forces 
which keep the condensate elongated even during the expansion, as shown for example 
by the red curve in figure 18(c). 

6.3. Collapse dynamics 

The dynamics of a condensate with pure contact interactions when the scattering 
length is ramped to a negative value, thus making the condensate unstable, is 
extremely rich: one observes a fast implosion ("collapse") of the condensate, followed 
by inelastic losses and a subsequent 'explosion' of the remnant condensate accompanied 
by energetic bursts of atoms. The behaviour of those "Bose-Novae" has been 
extensively studied with *^Rb and ^Li condensates [104, 131, 132, 133]. More recently, 
the formation of soliton trains during collapse has been reported [134, 135]. It is then 
natural to ask whether the collapse of a dipolar condensate displays some specific 
features arising from the long range and anisotropic character of the dipole-dipole 
interaction. 

In [136], the collapse dynamics of a dipolar ^^Cr BEG when the scattering length 
a is decreased (by means of a Feshbach resonance) below the critical value for stability 
Ocrit was investigated experimentally. A BEG of typically 20,000 atoms was created 
in a trap with frequencies {i^x, i^y, i^z) — (660, 400, 530) Hz at a magnetic field ~ 10 G 
above the 589 G Feshbach resonance, where the scattering length is a ~ 0.9 Obg- The 
scattering length a was then ramped down rapidly to a value af = 5ao, which is 
below the collapse threshold acrit — 13ao. After the ramp, the system evolved for 
an adjustable time ihoid and then the trap was switched off. The cloud was then 
imaged after time of ffight. The atomic cloud had a clear bimodal structure, with a 
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Figure 19. (a) Experimental images of a dipolar condensate after collapse and 
explosion, as a function of the time thold between the crossing of the critical 
scattering length for instability and the release from the trap. The time of flight 
is 8 ms. (b) Results of a numerical simulation of the collapse dynamics, without 
any adjustable parameter. The field of view is 130 fxra x 130 /Lim. 



broad isotropic thermal cloud, well fitted by a Gaussian, and a much narrower, highly 
anisotropic central feature, interpreted as the remnant BEG. Figur(^ 19(a) shows the 
time evolution of the condensate when varying thoid- From an initial shape elongated 
along the magnetization direction z, the condensate rapidly develops a complicated 
structure with an expanding, torus-shaped part close to the z = plane. Interestingly, 
the angular symmetry of the condensate at some specific times (e.g. at thoid = 0.5 ms) 
is reminiscent of the d— wave angiilar symmetry 1 — 3 cos^ 9 of the dipole-dipole 
interaction. Figure 19(b) displays the column density / |^(r)|^ Ax obtained from 
a numerical simulation of the three-dimensional GPE 



where 



2^A + Hrap + / t/(r-r',i)|^(r',t)|^dV- ^l^r 



V', (6.7) 



U{r,t) = ^-^^S{r) + mom! l^lScos^ 
m 47r r'^ 

stands for both contact and dipolar interactions. The non-unitary term proportional 
to L3 ~ 2 X 10"'*" m^/s describes three-body losses that occur close to the Feshbach 
resonance. The agreement between the experimental data and the simulation, 
performed without any adjustable parameter, is excellent. 

The observed cloverleaf patterns are caused by the anisotropic collapse and the 
subsequent dynamics of the system: when the atomic density grows due to the 
attractive interaction, three-body losses predominantly occur in the high-density 
region. The centripetal force is then decreased, and the atoms that gathered in 
this narrow central region are ejected due to the quantum pressure arising from the 
uncertainty principle. The kinetic energy is supplied by the loss of the negative 
interaction energy. As the collapse occurs mainly in the x — y direction due to 
anisotropy of the dipole-dipole interaction (in the absence of inelastic losses, the 
condensate would indeed bec;ome an infinitely thin cigar-shaped cloud along z, see 
section 5.3), and therefore the condensate "explodes" essentially radially, producing 
the anisotropic shape of the cloud. 

During the collapse, the BEG atom number, which was initially iVBEc(O) — 
16,000, dropped to a value ~ 6,000. The missing atoms very likely escaped from 
the trap as energetic molecules and atoms produced in three-body collisions. The 
simulated atom number as a function of thoid matched very well the experimental 
data. 



CONTENTS 



35 



The numerical simulation gives access not only to the density [-(/^(r)!^, but also 
to the phase S{r) of the order parameter tp {i.e. to the velocity field v = KVS/m) 
and reveals the generation of vortex rings. Figure 20(a) shows the simulated in- 
trap iso-density surface of a condensate at ihoid = 0.8 ms and the location of the 
vortex rings (shown as red curves). The mechanism responsible for the formation of 
vortex rings can be understood intuitively as follows. During the collapse, due to the 
strong anisotropy of the dipole-dipole interaction, the atoms ejected in the x — y plane 
flow outward, while the atoms near the z axis still flow inward, giving rise to the 
circulation. Thus, the vortex-ring formation is specific to the d-wave collapse induced 
by the dipole-dipole interaction. Although the vortex rings are not observed directly 
in the experiment (even when reconstructing the 3D density distribution by means of 
the inverse Abel transform), the excellent agreement between the experiment and the 
simulations strongly suggests the creation of vortex rings during the collapse. 

In [137], the collapse dynamics of a dipolar BEC was studied for different trap 
geometries, form prolate to oblate traps. The latter being created by superimposing a 
large period optical lattice onto the optical trap, it was possible to prepare independent 
condensates, let them collapse by changing the scattering length, and then release the 
confinement. The observation of high contrast interference fringes after the clouds 
overlapped (see figure 21) proved for the first time that the post-collapse remnant 
clouds are truly coherent matterwave fields. 

The collapse observed here has the same physical origin as the phonon instability 
discussed in section 5.1, and should be distinguished from other possible collapse 
mechanisms, in particular the roton instability discussed in section 5.6. Note also 
that the phonon instability in two-dimensional geometries [138] does not necessarily 
lead to collapse, as we discuss in section 7.3. 

7. Non linear atom optics w^ith dipolar gases 

As mentioned in section 4, BEC physics is inherently nonlinear due to the interparticle 
interactions. Non-dipolar BECs obey the non-linear Schrodinger equation (NLSE) 
(4.2), which is identical to that appearing in nonlinear optics of Kerr media. Striking 




Figure 20. Iso-density surfaces of a dipolar BEC after collapse and explosion, 
(a), simulation result, showing the location of vortex rings (in red), for the 
conditions of [136]. (b), experimentally reconstructed iso-density surfaces 
(using the inverse Abel transform) for high and low density (top and bottom, 
respectively) in the case of a cylindrically symmetric situation [137]. 




Figure 21. Interference pattern of independent condensates for different holding 
times thoid [137]. 

resemblances between both fields have been observed, including nonlinear atom- 
optics phenomena as four- wave mixing [139], BEC collapse [104], and the creation 
of bright [134, 140], dark [141, 142, 143] and gap [144] solitons. Dipolar BECs 
obey the NLSE (4.3), where the nonlinearity is intrinsically nonlocal, due to the 
long-range character of the dipole-dipole interaction. Nonlocality appears in many 
different physical systems, including plasmas [145], where the nonlocal response is 
induced by heating and ionization, nematic liquid crystals, where it is the result of 
long-range molecular interactions [146], and also photorefractive media [147]. Most 
interestingly, nonlocality plays a crucial role in the physics of solitons and modulational 
instability [148, 149]. In this section, we review some recent results concerning the 
nonlinear atom optics with dipolar BECs, including qualitatively new phenomena in 
bright and dark solitons (for recent works on vector and discrete solitons see [150] 
and [151]), vortices and pattern formation. 

7.1. Solitons 

Since the seminal works of Zakharov [152] it is known that the ID NLSE with 
focusing local cubic nonlinearity supports the existence of localized waves that travel 
with neither attenuation nor change of shape due to the compensation between 
dispersion and nonlinearity. These so-called bright solitons occur in diverse fields, 
most prominently in nonlinear optics [153]. Matter-wave bright solitons have been 
observed in quasi-lD condensates with a < [134, 140]. The quasi-lD condition 
requires a tight transversal harmonic trap of frequency uj± such that hu!± exceeds the 
mean-field interaction energy. This in turn demands the transversal BEC size to be 
smaller than the soliton width. When this condition is violated the soliton becomes 
unstable against transversal modulations, and hence multi-dimensional solitons are 
not stable in non-dipolar BECs. 

Remarkably the latter is not necessarily true in the presence of nonlocal 
nonlinearity. In particular, any symmetric nonlocal nonlinear response with positive 
definite Fourier spectrum has been mathematically shown to arrest collapse in 
arbitrary dimensions [149]. Multidimensional solitons have been experimentally 
observed in nematic liquid crystals [154] and in photorefractive screening media [155], 
as well as in classical ferrofluids [156]. Multidimensional solitons have been also 
discussed in BECs with short-range interactions, by considering the collapse inhibition 
induced by the first nonlocal correction to the local pseudopotential [148, 157]. 
However, this occurs for an extremely small BEC size [157], which, except for the 
case of a very small particle number, leads to extremely large densities, at which 
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Figure 22. Breathing (bold solid) and m = ±2 quadrupolc (bold dashed) mode 
of a 2D bright soliton for g = 20, with ko = Ip/lz the ratio between the condensate 
widths perpendicular to and along the dipole orientation. Results from a reduced 
2D NLSE [158] are shown in thin lines. Inset: E{kq) for g = 500, and = —0.42 
(dashed) and eja = —0.84 (solid). 

three-body losses destroy the BEC. 

Dipolar BECs on the contrary introduce a nonlocality at a much larger length 
scale. As a consequence, and in spite of the fact that, due to the anisotropy of the 
dipole-dipole interaction, the non local nonlinear response is not positive definite, 2D 
bright solitons may become stable under appropriate conditions [158]. This may be 
understood from a simplified discussion where we consider no trapping in the xy- 
plane and a strong harmonic confinement with frequency in the ^;-direction, along 
which the dipolcs arc oriented. A good insight on the stability of 2D solitons may be 
obtained from a Gaussian ansatz \E'(r) oc exp(— p-^/2?^ — z^/2Z^), where Iz = \Jh/mijJz, 
and Ip = kqIz is the a;y-width. Introducing this ansatz into the nonlocal NLSE (4.3) 
we obtain the system energy, which up to a constant is 

i;=^{27r + 5[l-edd/(«o)]} (7.1) 

where g = g/\/2TTlTLJzll, and /(k) is defined in (5.10). As mentioned above, in the 
absence of dipole-dipole interaction (sdd = 0), 2D localized solutions are unstable, 
since E{k,q) oc Kq"^ cither grows with Lp (collapse instability) or decreases with Lp 
(expansion instability). Dipolar BECs are remarkably different due to the additional 
dependence /{kq), which may allow for a minimum in E{ko) (inset in figure 22), i.e. for 
a stable localized wavepacket. From the asymptotic values /(O) = —1 and /(oo) = 2 
localization is just possible if edd5 < 27r + g < — 2edd.9- A simple inspection shows that 
this condition is fulfilled only if Sdd < 0, i.e. for the dipole-dipole interaction tuned 
with rotating fields (see section 2.2). Note also that if Na/l^ ^ 1, the stability 
condition reduces to |edd| > 1/2. The anisotropic character of the dipole-dipole 
interaction becomes particularly relevant when relaxing the quasi-2D condition. In 
particular, solitons in 3D dipolar BECs are fundamentally unstable against collapse, 
as reflected in the decrease of the frequency of the breathing mode of the soliton for 
larger values of £dd [158] (figure 22). 

Although Sdd < is attainable for magnetic dipolcs in a rotating magnetic 
field, the combination with Feshbach resonances to reduce the contact interactions 
makes experiments, e.g. in Chromium, very complicated. Recently a similar idea, but 
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Figure 23. Stability diagram of an anisotropic soliton as a function of and 
9cr = Qci/V^lz, where for g > ga the soliton is unstable against collapse even 
for Sdd > 1- 




Figure 24. Density plot of the fusion of two dipolar 2D solitons for g = 20, 
Sdd = ~2.1, and initial relative momentum kl^ = 0.01 along x. From top to 
bottom ,uj^t/2000 = 0, 1, 2, 3, 4, and 5. 



without the necessity of dipole tuning, has been proposed by Tikhonenkov, Malomed 
and Vardi [159]. In this new proposal the dipoles are considered as polarized on the 
2D plane. A similar Gaussian ansatz on the xy-plane as above with unequal widths 
Ix and ly shows the appearance of a minimum in the energy functional E{lx,ly) for 
£dd > 1; and thus the existence of stable anisotropic solitons. However, even if the 
latter condition is fulfilled there is a critical universal value ^cr (which decreases with 
£dd) such that for g > ^cr the minimum of E[lx, ly) disappears [160] (see figure 23). As 
a consequence, contrary to the case of isotropic solitons (with Sdd < 0) there is a critical 
number of particles per soliton even if the soliton remains 2D. In addition, as for the 
case of isotropic solitons, anisotropic solitons are also unstable in 3D environments. 

A major difference between bright solitons in non-dipolar and dipolar BECs 
concerns the soliton-soliton scattering properties. Whereas solitons in ID non-dipolar 



CONTENTS 



39 



BECs scatter elastically, the scattering of dipolar solitons is inelastic due to the 
lack of intcgrability [148]. The solitons may transfer ccntcr-of-mass energy into 
internal vibrational modes, resulting in intriguing scattering properties, including 
soliton fusion [158] (see figure 24), the appearance of strong inelastic resonances [161], 
and the possibility of observing 2D-soliton spiraling as that already observed in photo- 
refractive materials [155]. 

For defocusing nonlinearity (a > 0) the local NLSE supports dark-soliton 
solutions, i.e. density notches (accompanied by phase slips) that propagate 
with no change of shape, again due to the compensation between dispersion 
and nonlinearity [162]. Dark solitons have been created in non-dipolar quasi-lD 
BECs [141, 142, 143], but become fimdamcntally imstablc in higher dimensions against 
vibrations of the nodal plane which lead to the so-called snake instability. This 
instability which was previously studied in the context of nonlinear optics [163, 164] 
leads in the context of non-dipolar BEC to the soliton break-down into vortex rings 
and sound excitations [165, 166, 167, 168]. On the contrary, a dark soliton in a 
dipolar BEC may become stable in a 3D environment [138] if the BEC is placed in a 
sufficiently deep 2D optical lattice (characterized by an effective mass on the lattice 
plane m* > m [169] and a regularized local coupling constant g). This effect may 
be understood by considering the phonon-like excitations of the dark soliton plane: 
e = \J a jMq^ where M is the (negative) soliton mass per unit area, and ct is the 
surface tension of the nodal plane. For non-dipolar BECs one always has a > 0, 
and hence the phonon spectrum is always unstable leading to the above mentioned 
snake instability. However, for dipolar BECs the surface tension a becomes negative 
when m/m* < 3Cdd/(35 -|- 2Cdd) [138], and hence the nodal plane is stabilized for 
a sufficiently deep optical lattice (large m*/m) and a sufficiently large dipole-dipole 
interaction. 

7.2. Vortices 

Quantized vortices constitute one of the most important consequences of superfluidity, 
playing a fundamental role in various physical systems, such as superconductors [170] 
and superfluid helium [171]. Vortices and even vortex lattices have been created in 
BECs in a series of milestone experiments [172, 173, 174]. Contrary to normal fluids, 
a vortex in a BEC cannot be created by any rotation, but there is a critical angular 
velocity Qc- Only for angular velocities f2 > fie it is energetically favorable to form 
a vortex [175] (the critical rotation for vortex nucleation is actually larger, since also 
dynamical instabilities at the condensate boundaries must be considered [176]). O'Dell 
and Eberlein [177] have recently studied the critical f2c in a dipolar condensate in the 
Thomas-Fermi regime, by means of the solution discussed in section 5. For BECs in 
axially symmetric traps with the axis along the dipole orientation it has been shown 
that Qc is decreased due to the dipole-dipole interaction in oblate traps and increased 
in prolate traps. This modification can be traced back to the modification of the 
Thomas- Fermi radius due to the dipole-dipole interaction, rather than changes in the 
vortex core. However, the dipole-dipole interaction may induce crater-like structures 
close to the vortex core for the case of a < 0, or even anisotropic vortex cores, as 
recently shown by Yi and Pu [178]. 

At higher rotating frequencies, more vortices enter the condensate and a vortex 
lattice develops [173, 174]. In non-dipolar BECs vortices form so-called Abrikosov 
lattices, i.e. triangular lattices with hexagonal symmetry. Interestingly, this is not 
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Figure 25. Stable/unstable regimes, as a function of the ratio m/m* and e^d, 
for straight vortex lines when the dipoles are oriented along the vortex line. 



necessarily the case in dipolar condensates [179, 180]. In particular, with increasing 
Edd or high filling factor, the vortex lattice may undergo transitions between different 
symmetries: triangular, square, stripe vortex crystal, and bubble states [179]. In 
addition for vortex lattices in double well potentials the competition between tunneling 
and interlayer dipole-dipole interaction should lead to a quantum phase transtion from 
a coincident phase to a staggered one [180]. 

Vortex lines are in fact 3D structures with transverse helicoidal excitations known 
as Kelvin modes [181, 182]. Kelvin modes, which play an important role in the physics 
of superfluid Helium [171], have also been experimentally observed in BECs [183]. 
The long-range character of the dipole-dipole interaction links different parts of the 
vortex line, and hence the 3D character of the vortex lines is much more relevant in 
dipolar BECs. Remarkably the dipole-dipole interaction may significantly modify the 
vortex-line stability. In the presence of an additional optical lattice (which leads to 
an effective mass m* along the vortex line direction) the dispersion of Kelvin modes 
shows a roton-like minimum [184], which for sufficiently large dipole-dipole interaction 
and large m* may reach zero energy, leading to a thermodynamical instability related 
to a second-order-like phase transition from a straight vortex into a twisted vortex 
line [185] (see figure 25). 

7. 3. Pattern formation 

Pattern formation in driven systems is a general nonlinear phenomenon occurring in 
many scenarios ranging from hydrodynamics and nonlinear optics to liquid crystals and 
chemical reactions [186]. Faraday patterns have been recently observed in non-dipolar 
BECs by a modulation of the harmonic confinement [187] which leads to a periodic 
modulation of the system nonlinearity [188]. Faraday patterns offer an important 
insight about the elementary excitations in a BEC, and hence pattern formation may 
be remarkably different in dipolar BECs, especially in the presence of a roton-maxon 
excitation spectrum (section 5). Most remarkably, whereas for non-dipolar BECs the 
pattern size decreases monotonously with the driving frequency, patterns in dipolar 
BECs present a highly non-trivial dependence characterized by abrupt pattern-size 
transitions [115]. 

Faraday pattern formation in driven systems is the (transient) result of an 
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externally driven dynamical instability. However, an interaction-induced dynamical 
instability may lead as well to pattern formation. In this sense, the phonon instability, 
which, as discussed in section 5, leads to collapse in 3D dipolar gases (also in 2D and 
3D non-dipolar BECs), does not necessarily lead to collapse in 2D geometries. On 
the contrary, the 2D phonon instability leads to the formation of a soliton gas, and to 
transient pattern formation, which, if avoiding collapse, may lead to the formation 
of a 2D stable soliton [138]. Dynamical roton instability leads typically to local 
collapses [111], although a sufficiently strong trapping may stabilize a biconcave BEG 
profile [100] (as mentioned in section 5.3). 

8. Dipolar effects in spinor condensates 

Spinor BECs, composed of atoms in more than one Zeeman state, constitute an 
extraordinary tool for the analysis of multicomponent superfluids. Whereas magnetic 
trapping confines the BEG to weak-field seeking magnetic states, optical trapping 
enables confinement of all magnetic substates, hence freeing the spin degree of 
freedom [189]. Interestingly, interatomic interactions allow for a coherent transfer 
of population between different Zeeman states (spin-changing collisions), leading 
to a fascinating physics in both what concerns ground-state properties and spin 
dynamics [190, 191, 192, 193, 194, 195]. 

The energy scale associated to spin-preserving collisions is given by the chemical 
potential, which for typical alkali gases (and even for Ghromium in absence of Feshbach 
resonances) is much larger than the dipole-dipole interaction. On the contrary the 
energy associated with spin-changing collisions is typically much smaller, since it 
is provided by the difference between s-wave scattering lengths in different spin 
channels [190], which is very small. Hence the dipole-dipole interaction may become 
comparable to the energy of spin-changing collisions, and as a consequence even alkali 
spinor BEGs (in particular ^''Rb) can be considered in this sense as dipolar gases as 
well. The dipole-dipole interaction may hence play a significant role in the ground- 
state properties and the dynamics of spinor condensates. 

8.1. Ground state 

Pu et al. [196] have shown that ferromagnetic spinor BEGs (as it is the case of ^^Rb 
in f = 1) placed at different sites of a strong optical lattice behave as single magnets 
oriented in the effective magnetic field induced by the combination of a external 
magnetic field and the dipole-dipole interaction of other sites. Interestingly, the 
collectively enhanced magnetic moment of the condensates at each site enhances the 
magnetic dipole-dipole interaction between sites, which may become sufficiently strong 
even for alkali atoms [196]. As a consequence, such an array of effective magnets 
can undergo a ferromagnetic (ID lattice) or anti- ferromagnetic (2D lattice) phase 
transition under the magnetic dipolar interaction when external magnetic fields are 
sufficiently weak [196, 197]. In addition, for ID lattices the inter-site dipole-dipole 
interaction may distort the ground-state spin-orientations and lead to the excitation 
of spin waves, which possess a particular dispersion relation which depends on the 
transverse width of the condensates [198]. 

The dipole-dipole interaction may play an important role in the properties of 
trapped spinor BEGs, especially in the absence of significant external magnetic fields. 
In particular, whereas in absence of an external magnetic field the spinor BEG is 
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rotationally invariant in spin space, the dipole-dipole interaction breaks this symmetry, 
inducing new quantum phases which can be be reached by tuning the effective strength 
of the dipole-dipole interaction via a modification of the trapping geometry [199]. For 
the case of spin-1 BECs, for very low magnetic fields (typically below 10 jiG) the 
phase diagram presents due to the dipole-dipole interaction three different ground 
state phases, characterized by different spin textures: a polar-core vortex phase, a 
flower phase and a chiral spin- vortex phase [200]. The latter has chirality in the 
formation of the spin vortex, and the topological spin structure spontaneously yields 
a substantial net orbital angular momentum. 

8.2. Dynamics and Einstein- de Haas effect 

The dipole-dipole interaction has been shown to play quite a small role in the ground 
state properties of Chromium BECs [201, 202]. However the effects of the dipole- 
dipole interaction on the spinor dynamics may be much more intriguing. The short- 
range interactions in a spinor condensate may occur in different scattering channels, 
corresponding to different total spin of the colliding pair [190] 

2F 

Vsr = \ UrY^gsMr), (8.1) 

where Vs is the projector on the total spin S (necessarily even due to symmetry 
reasons), gs = ATTfi^as/m, and as is the s-wave scattering length for the channel of 
total spin S. The short range interactions necessarily preserve the spin projection Sz 
along the quantization axis. 

The dipole-dipole interaction for a spinor BEC is of the form 

- 3{Smn ■ e){Sm'n' ' e)\ iJn{r)iJn' {v'), (8.2) 

where S = {Sx,Sy,Sz) are the spin-F matrices, and Cdd = /ioA*Bfi'|'/47r (for ^^Cr 
F = 3 and dd = 0.004^6), with e = (r - r')/\r - r'\. 

Interestingly, contrary to the short-range interaction the dipole-dipole interaction 
docs not necessarily conserve the spin projection along the quantization axis due to 
the anisotropic character of the interaction. In particular, if the atoms are initially 
prepared into a maximally stretched state, say mp = —F, short-range interactions 
cannot induce any spinor dynamics due to the above mentioned conservation of 
total magnetization S^. Dipole-dipole interactions, on the contrary may induce a 
transfer into mp + If the system preserves cylindrical symmetry around the 
quantization axis, this violation of the spin projection is accompanied by a transfer of 
angular momentum to the center of mass, resembling the well known Einstein-de Haas 
effect [202, 203]. Due to this transfer an initially spin-polarized dipolar condensate 
can generate dynamically vorticity (see figure 26). 

Unfortunately, the Einstein-de Haas effect is destroyed in the presence of even 
rather weak magnetic fields. Typically, magnetic fields well below 1 mG are necessary 
to observe the effect. Due to the dominant role of Larmor precession, and invoking 
rotating- wave-approximation arguments, the physics must be constrained to manifolds 
of preserved magnetization, where the system presents a regularized dipole-dipole 
interaction [204]. However the dipole-dipole interaction may have observable effects 
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Figure 26. Isodensity surfaces of mp = +1, nip = 0, and mp = — 1 spinor 
components corresponding to the densities: 7.96 X 10^"^ cm~^, 7.96 X 10^^ cm~'^, 
and 7.96 X 10^^ cm~'^, respectively. The color on the surface represents the phase 
of the component wave function (with the scale given on the right). Note the 
characteristic vortex rings patterns. The trap frequency uJx,y,z = 2tt X 100 Hz and 
the magnetic field B = — 0.029 mG. The density plots are taken at 140 ms. 

also under conserved magnetization, since the regularized dipole-dipole interaction 
may lead also in that case to spin textures [204] . 

A significant Einstein-de Haas effect may be recovered however under particular 
resonant conditions as studied by Gawryluk et al. [205], who studied the particular 
case of Rb BECs in _F = 1, initially prepared into mp = 1. In that case, the population 
transfer away from the mp — I state is typically very small, but it can be significantly 
enhanced by applying a resonant magnetic field, such that the Zeeman energy of an 
atom in the mp = ^ state is totally transferred into the kinetic energy of the rotating 
atom in tof = 0, fiB — i?kin- Typically, the resonant B is small (~ 100/iG), but 
can be tuned directly, or by adjusting the trap geometry and thus the rotational 
energy of the atoms, reaching up to 1 mG. Gawryluk et al. demonstrated that at the 
resonance, a significant transfer of the initial population of the mF = 1 state occurs 
on a time scale inversely proportional to the dipolar energy ttransfcr — fi/C'ddn. As 
mentioned above, such a transfer is accompanied by the formation of vorticity (see 
figure 26) [202, 203, 205]. 

An alternative possibility for the observation of the Einstein-de Haas effects at 
finite fields may be also provided by an artificial quadratic Zeeman effect induced by 
either microwave [206] or optical fields [207]. This effective quadratic Zeeman effect 
may allow for a resonance between e.g. mp — —F and mp — ~F + 1 and hence lead 
to a significant enhancement of the Einstein-de Haas transfer [207] . 

8.3. Experimental results 

Recent 2D experiments at Berkeley [208] have demonstrated the dipolar character 
of spin-1 ^''Rb spinor BECs. In particular, these experiments show the spontaneous 
decay of helical spin textures (externally created by magnetic field gradients) toward 
a spatially modulated structure of spin domains (see figure 27) . The formation of this 
modulated phase has been ascribed to magnetic dipolar interactions that energetically 
favor short-wavelength domains over the long-wavelength spin helix. Interestingly, 
the reduction of dipolar interactions (by means of radio-frequency pulses) results in a 
suppression of the modulation. 

These experiments have attracted a large deal of theoretical interest, in particular 
in the dipole-dipole interaction-induced distortion of the excitation spectrum of a 
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Figure 27. Spontaneous dissolution of helical textures in a quantum degenerate 
^^Rb spinor Bose gas. A transient magnetic field gradient is used to prepare 
transversely magnetized (b) uniform or (a),(c) helical magnetization textures. The 
transverse magnetization column density after a variable time T of free evolution is 
shown in the imaged x — z plane, with orientation indicated by hue and amplitude 
by brightness (color wheel shown), (b) A uniform texture remains homogeneous 
for long evolution times, while (c) a helical texture with pitch A = 60 fim dissolves 
over ~ 200 ms, evolving into a sharply spatially modulated texture. Figure 
reprinted with permission from [208]. Copyright 2008 by the American Physical 
Society. 

spinor BEC [209, 210]. Recently Cherng and Demler have analyzed (in the context of 
the above mentioned 2D experiments [208] ) the possibility of roton softening (similar 
to that discussed in section 5.6) in the spectrmii of spin excitations. This roton 
instability may lead as a function of the quadratic Zeeman effect and the magnetic 
field orientation with respect to the normal of the 2D trap to different spin textures 
(checkerboard, striped phase) and may pave the way towards a supersolid phase. 
The latter is a long pursued phase in condensed-matter physics which possesses both 
periodic crystalline order and superfluidity [211, 212, 213] (see also the next section). 
However, instability against finite-momentum excitations does not necessary lead to 
the appearance of a stable modulation (as it is the case of the roton instability 
discussed in section 5.6) and more work on the physics in the unstable regime is clearly 
necessary. Recent experiments at Berkeley have reported on possible first traces of 
supersohdity in a spinor ^''Rb BEC [214]. 

9. Dipolar gases in optical lattices 

One of the most fruitful fields of research both from the experimental and the 
theoretical points of view in the last years has been the study of ultra-cold atomic 
samples in optical lattices, which are non dissipative periodic potential energy surfaces 
for the atoms, created by the interference of laser fields. The study of cold atoms 
in periodic potentials is of primary interest because it allows to reproduce problems 
traditionally encountered in condensed matter and solid state physics in a new setting, 
where a high degree of control is possible and where the Hamiltonian which governs the 
system is in general very close to some idealized one. With the present developments 
there is even the possibility to investigate, with appropriately designed atomic systems, 
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phenomena which do not exist in condensed matter. 

We will first summarize the physics of weakly interacting atomic gases in optical 
lattices in section 9.1 and describe the first measurement of dipolar effects in alkali 
atoms in section 9.2. After introducing the physics of strongly correlated systems 
with contact interaction in section 9.3, we will devote section 9.4 and section 9.5 
respectively to the quantum phases and the metastable states found in 2D lattices in 
the presence of long-range interactions. In section 9.6, we discuss the novel physics 
introduced by the presence of two or many 2D optical lattice layers. In section 9.7, 
we discuss proposals on how to tailor the interaction potential and create lattice spin 
models with polar molecules and finally, in section 9.8, we will mention the possibility 
of formation of self- assembled regular structures in cold dipolar gases. 

9.1. Bose-Einstein condensates in optical lattices 

In the presence of weak optical lattices, when the coherence of the system is preserved, 
the Gross-Pitaevskii equation provides a good description of the system. Due to the 
presence of the periodic potential and interactions, analogies to phenomena typical of 
solid state physics and non-linear optics are respectively possible [215, 216, 217, 14]. 

As previously explained, in the most common cases, interactions in ultra-cold 
gases are dominated by s-wave scattering, which can be in a very good approximation 
considered a point-like interaction. In the case of full coherence, the system is described 
by a macroscopic wavefunction and obeys the GPE (4.5). In the presence; of optical 
lattices, T4xt is periodic and given by I4xt(^) = S„ Ki'^* sin^(7ra;„/d„), where the 
index n runs over the dimensions of the lattice and dn is the lattice constant in the 
n-th direction. For lattices created by counterpropagating laser beams of wavelength 
A, the lattice spacing is d = A/2. The depth of the lattice potential V"^^ depends 
linearly on the intensity of the laser light. 

It is well known that the spectrum of a single particle in a periodic potential is 
characterized by bands of allowed energies and energy gaps [218]. The counterpart 
of the energy eigenstates delocalized over the whole lattice (Bloch states) are the 
wavefunctions centered at the different lattice sites (Wannier functions). For deep 
enough periodic potentials, when the Wannier functions are well localized at the 
lattice sites, the so-called tight binding regime is reached: the first band takes the 
form E{q) — —2 cos((7„c?„), g„ being the quasi-momenta in the different lattice 
directions and J the tunneling parameter between neighboring wells. For low enough 
interactions and temperature, the physics of the system is well approximated by the 
one taking place in the first energy band. Under those assumptions, the discretized 
non-linear Schrodinger (DNLS) equation [219, 220, 221, 222, 223] provides a good 
description of the system. 

The excitations on top of the GP solution can be found by generalizing the 
Bogoliubov method to include the periodic potential. They show a phononic branch 
with renormalized sound velocity [224, 225, 226]. For a moving condensate, one finds 
complex frequencies at the edge of the Brillouin zone [227] , highlighting the presence 
of dynamical instabilities [228]. 

Among the numerous phenomena described theoretically and observed 
experimentally are collective oscillations [229], Bloch oscillations [230], dynamical 
instabilities [231], Josephson oscillations [232], non linear self-trapping [233], gap 
solitons [144]. The collective oscillations can be described through an effective 
macroscopic dynamics accounting for the periodic potential and their frequency is 
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rescaled in terms of the effective mass [169]. Instead the physics underlying the other 
phenomena is dominated by a crucial interplay between the periodic potential and 
interactions. 

In the alkalies usually used in experiments with optical lattices, s-wave scattering 
dominates over all other typos of interactions. When the s-wavc scattering length 
is reduced, e.g. by means of a Feshbach resonance, the presence of other types of 
interaction is relatively enhanced and can be probed. In the next section, we present 
the first measured effects of dipole-dipolc interaction for a quantum gas in an optical 
lattice, performed with Potassium atoms. 

9.2. Block oscillation damping due to dipole-dipolc interactions 

Bloch oscillations have been one of the first solid state phenomena to be investigated 
with cold atoms [234, 235], observed shortly after Bloch oscillations of electrons in 
semiconducting heterostructures [236, 237, 238]. They consist in oscillations in space 
in the presence of a constant acceleration, due to the change of sign of the effective 
mass along the first energy band. For electrons, the acceleration is provided by a 
constant electric field, while for cold atoms it is produced by a linear increase in time 
of the relative detuning of the two laser beams creating the optical lattice or, in a 
vertical setup, by gravity. 

With Bose-Einstein condensates, Bloch oscillations can be measured with a higher 
precision thanks to the smaller width of the momentum distribution [230]. However 
in the presence of interatomic interactions, the on-set of dynamical instabilities in the 
outer region of the Brillouin zone causes a damping of the oscillations. Due to the 
relevance of Bloch oscillations as a tool of precision measurement of accelerations and 
small forces at small distance from surfaces [239, 240, 241, 242], attempts of reducing 
interactions have been pursed first with polarized fermions [243] and then with bosons 
with reduced s-wave scattering length [244, 128]. 

The experiment in [128] has highlighted for the first time effects of dipole- 
dipole interactions on alkali atoms and has shown how its interplay with the contact 
interaction can be exploited to reduce interaction-induced decoherence of Bloch 
oscillations in a ID optical lattice. When the scattering length is tuned to zero by 
means of a Feshbach resonance, the magnetic dipole-dipole interaction becomes the 
limiting factor for the coherence time. The point of minimum decoherence is shifted 
to negative or positive values of the scattering length, depending on the orientation 
of the dipoles with respect to the axis of the lattice (see figure 28). 

The values of the scattering lengths which maximize the lifetime in the different 
configurations have been predicted through the solution of ID DNLS and GP 
equations including the dipolar potential. The 3D geometry has been accounted for, 
assuming the condensate to be in the transversal groimd state. The dipole-dipole 
interaction contributes both to the regularization of the on-site interactions and also 
to the intersite interactions. Although the compensation of the on-site dipole-dipole 
interaction and the short-range interactions explains the sign of the displacement 
(towards a < or o > 0) of the decoherence minimum, the dipole-induced intersite 
interactions are crucial for the quantitative understanding of the experimental results 
in [128]. In this sense the experiments in [128] constitute also the first observation of 
intersite effects in dipolar gases in optical lattices. Due to the intersite interactions, 
the dipole-dipole interaction can be only partially compensated by the short-range 
interactions, and as consequence a non complete reduction of the decoherence rate of 
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Figure 28. Decoherence of Bloch oscillations in "^^K around the scattering 
length zero crossing, due to the interplay of contact and dipolar interactions. 
The character of the dipolar interaction depends on the relative orientation of the 
lattice and magnetic field; (a) prevalently repulsive interaction; (b) prevalently 
attractive interaction, (c) The width of the momentum distribution after a few 
100 ms of Bloch oscillations shows a minimum when the contact interaction 
compensates the dipolar one: circles are for the case (a), while squares are for the 
case (b). The lines are parabolic fits to the data, which constrain the position of 
the zero-crossing in a comparison with theory (black region in the lower panel) 
better than Feshbach spectroscopy (gray region). Figure courtesy of G. Modugno. 

the ^^K-based interferometer is observed (with a minimum residual rate of 0.05 Hz). 

In the experiments in [128], the interactions at the decoherence minimum were 
much weaker than the transversal confinement, justifying the approximation of 
assuming the BEC in the transversal ground state. This is however not the general 
case. For shallower transversal confinements the intersite dipole-dipole interaction 
may significantly modify and even destabilize the spectrum of elementary excitations 
of a BEC in an optical lattice [114, 112]. In particular, the intersite interactions may 
induce rotonization and even roton-instability (see also section 5.6) under appropriated 
conditions, and may (for a sufficiently shallow transversal confinement) lead to a 
dynamical instability that could as well damp Bloch oscillations [114]. Interestingly, 
the intersite dipole-dipole interaction induces an hybridization of transversal modes 
at different sites (and a corresponding band-like spectrum) even if the hopping is 
completely suppressed [112]. Remarkably, whereas a single lattice site could be 
stable, a stack of non-overlapping dipolar BECs may become roton-unstable, showing 
once more that polar gases in optical lattices differ qualitatively from short-range 
interacting gases. 

9.3. Strongly correlated lattice gases. Bose-Hubbard Hamiltonian 

For deep optical lattices and small numbers of atoms per site, the coherent description 
of the system provided by the GPE breaks down due to the growing effect of 
correlations. One of the greatest achievements of the last years was the experimental 
observation [16] of the superfluid to Mott insulator transition [245, 15]. In this session, 
we will introduce the main theory for ultra-cold atoms in optical lattices in the case of 
point-like interaction, providing the background for the case of long-range interactions. 
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which will be treated in the following sections. 

Under usual experimental conditions, the single band approximation mentioned 
in section 9.1 is appropriate. In order to allow for the breaking of the coherence of 
the system, the field operator is replaced by its single-band many-mode expansion 
'4'{r) = ^jWi(r)ai, with cii being the annihilation operator for one boson in the 
Wannier function Wi{r) localized at the bottom of the lattice site i. 

Neglecting the overlap beyond nearest neighboring densities, defining 

J =- j <ir) (^-^ + yext(r)) Wi+^{r) dV, (9.1) 

U = 9 j ki(r)|^ dV, (9.2) 
and rij = ajaj, one can derive [15] the famous Bose-Hubbard Hamiltonian 



—ni{ni - 1) - iirii 



(9.3) 



(ij) i 

extensively studied in condensed matter physics. In optical lattices, the Hamiltonian 

parameters can be accurately controlled by changing the light intensity: ramping it 
up increases the interaction term U due to a stronger localization of the wavefunctions 
at the bottom of the lattice wells, and at the same time exponentially decreases the 
tunneling J. 

When tunneling is suppressed compared to interactions, this Hamiltonian presents 
a quantum phase transition between a superfluid phase, characterized by large number 
fluctuations at each lattice site, and a Mott insulating phase where each lattice well is 
occupied by precisely an integer number of atoms without any number fluctuations. 
The nature of this phase transition and the qualitative phase diagram can be inferred 
based on very simple arguments [245] . 

At zero tunneling J = and commensurate filling (exactly an integer number 
n of atoms per well), the interaction energy is minimized by populating each lattice 
well with exactly n atoms. Energy considerations tell that the filling factor n is 
energetically favored in the range of chemical potential (n — 1)U < ji < nU . The 
state with precisely integer occupation at the lattice sites is called Mott insulating 
state. Since a particle-hole excitation at J = costs an energy /S.E = U equal to the 
interaction energy, the Mott state is the lowest energy state at commensurate filling. 
For a tunneling J different from zero the energy cost to create an excitation decreases 
thanks to the kinetic energy favoring particle hopping. However, for large interactions 
and small tunneling, the gain in kinetic energy (~ J) is not yet sufficient to overcome 
the cost in interaction energy (^ U), which leads to the existence of Mott insulating 
states also at finite tunneling. For large enough tunneling, instead, particle hopping 
becomes energetically favorable and the system becomes superfiuid. The regions in 
the J vs. /i phase diagram where the Mott insulating state is the ground state are 
called Mott lobes (see figure 29). For non commensurate filling, there are extra atoms 
free to hop from site to site at no energy cost, so that the phase of the system is 
always superfluid. The superfluid phase at non commensurate densities survives down 
to J = for fj,/U = [p] where the symbol [p] indicates the integer part of the density. 
Due to the finite energy cost required to add or remove one particle, the Mott phase 
is gapped and incompressible, while in the superfluid regions the gap vanishes and the 
system is compressible. 
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Figure 29. Schematic phase diagram for the ground state of the Bose-Hubbard 
Hamiltonian (9.3). 



In order to find the shape of the lobes at finite J sophisticated calculations are 
required. Apart from the mean- field approximation [245, 246, 247], which works only 
qualitatively in one dimension and works better and better in larger dimensions, there 
is no exact analytical method which allows to calculate the boundary of the lobes. 
Improvements are achieved by high order perturbative strong coupling expansions 
[248, 249] and exact numerical results are obtained using Quantum Monte Carlo 
techniques (see e.g. [250, 251, 252]). 

In the experiments, the phase transition has been identified by looking at the 
interference of the expanded cloud and the measurement of the gapped excitations 
in the Mott phase [16]. The SF-MI shell structure which arises in the presence of 
an external confinement [15, 253, 254] has been observed using spatially selective 
microwave transition and spin-changing collisions [254] or using clock shifts [255] , and 
the underlying ordering in the lattice in the Mott phase has been inferred from the 
measurement the periodic quantum correlations in the density fluctuations in the cloud 
after expansion [256, 257]. The interference of the expanding clond and spatial noise 
correlations measurements prove to be useful tools to identify also the exotic quantum 
phases expected to appear in the presence of dipole-dipole interaction. 

5.^. Quantum phases of dipolar lattice gases 

Dipole-dipole interactions add to the Bosc^-Hubbard model a new essential ingredient, 
given by long-range and anisotropic interactions. For a lattice of polarized dipoles, 
as sketched in figure 30, the extended Bose-Hubbard Hamiltonian in the presence of 
long-range interaction is 



where £ is the distance between the two optical lattice sites i and j. 

The extended Bose-Hubbard Hamiltonian (9.4) has been extensively studied. It 
has been predicted that in 2D lattices the presence of finite range interactions (where 
the sum over f. is generally cut-off at the nearest or next-nearest neighbor) gives 
rise to novel quantum phases, like the charge-density wave (checkerboard), namely 
an insulating phase with modulated density, and the supersolid phase, presenting 
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Figure 30. Schematic representation of a gas of polarized dipolcs in a ID optical 
lattice (a) and in a single 2D optical lattice layer (b). 

the coexistence of superfluidity and of a periodic spatial modulation of the density, 
different from the one of the lattice [258, 259, 260] (see also discussion in section 5.6). 

The presence of insulating phases at fractional filling factors can be inferred 
readily following the criteria explained in the previous section, which predict that 
(e.g.) a checkerboard ordering of the atoms (see figure 31(GS)) at J = is stable 
against particle-hole excitations in the range of chemical potential < /i < 4J7nn7 C^nn 
being the first nearest-neighbor interaction (and neglecting for the sake of simplicity 
all following ones). Analogous to the standard Mott insulating phase, these insulating 
phases at fractional filling factor exist in some low tunneling region of the fi vs. J 
phase diagram. 

In [261], it has been pointed out that ID lattice systems of spinless bosons 
interacting with long range interactions, possess a further insulating phase, which 
they call Haldane Bose Insulator (HI), presenting some analogies with the famous 
Haldane gapped phase in quantum spin-1 chains [262]. This is a gapped phase, which 
unlike the checkerboard phase does not break the translational symmetry of the lattice, 
but is characterized by an underlying hidden order, namely a non trivial ordering of 
the fluctuations which appear in alternating order separated by strings of equally 
populated sites of arbitrary length. 

The existence of the supersolid phase in solid helium has not yet been 
unambiguously proven experimentally: while on the one hand the interpretation of the 
first experimental results measuring a non-classical rotational inertia [263, 264, 265] 
remains controversial, microscopic calculations [266] indicate that disorder-based 
mechanisms, like the presence of superfluid dislocations, grain boundaries, and ridges, 
should be responsible for the more recent observations of supersolidity [267, 268] . Even 
if bulk supersolid remains in many respects more challenging than lattice supersolid, 
the question of the stability of the supersolid phase in the presence of the lattice 
is not trivial and has been only recently settled by exact Quantum Monte Carlo 
simulations. Checkerboard supersolid (at p w 1/2) is expected for dominant nearest- 
neighbor interaction, while star (at p « 1/4) and striped (at p « 1/2) supersolids* are 
predicted for non vanishing nearest-neighbor interaction [269, 270]. There have been 

* The star supersolid shows modulations of the density and order parameter such that in a 2 X 2 
elementary cell, one [the other] diagonal contains sites with different [equal] density and order 
parameter; the striped supersolid presents alternating horizontal or vertical striped of higher and 
lower density and order parameter. 
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several studies devoted to the stability of the supersolid phase versus phase separation 
[271, 272], which is identified by a negative compressibility. Agreement seems to be 
reached on the conclusion that the checkerboard supersolid is stabilized at p > 1/2 
by a finite on-site interaction and a strong enough nearest neighbor interaction [273], 
while it phase separates for nearest-neighbor interactions at densities smaller than 
1/2 (unless a strong enough nearest-neighbor hopping is introduced [274]). Instead 
the striped supersolid, obtained for large next nearest-neighbor interaction, exists for 
all doping away from p — \/2 both in the hard-core and soft-core cases [271, 272, 274]. 
Finally, at large next nearest-neighbor interaction, the star supersolid can always be 
obtained by doping a star solid at p = 1/4 with vacancies and by doping it with 
bosons in the case the p = 1/2 ground state is a striped crystal [274, 275]. The 
most important conclusion on which most papers agreefj is that no supersolid phase 
is found at commensurate density. Analogous results have been recently discussed for 
ID geometries [277, 278], the most important difference being the absence of phase 
separation. 

Providing an alternative setting where to look for the supersolid phase, cold 

atoms with long range interactions are particularly appealing [269, 270, 273, 279, 
280, 261]. Dipolar atoms and polar molecules are good candidates to creates 
such physical systems, bringing into play the extra feature of the anisotropy of 
the interaction. Dipolar gases have been first identified as possible candidates to 
provide a long-range interaction system in [269]. The on-site parameter U in (9.4) 
is given by two contributions: one is arising from the s-wave scattering Ug = 
Awfi^a/m J n^(r) d^r, and the second one is due to the on-site dipole-dipole interaction 

Udip = l/(27r)^ / C/dd(Q)^^(Q) d^q, UddiQ) and h{q) being the Fourier transforms of 
the dipole potential and density, respectively [97]. Because of the localization of the 
wavefunctions at the bottom of the optical lattice wells, the long range part of the 
dipole-dipole interaction Ue. is in a very good approximation given by the dipole-dipole 
interaction potential at distance ^, Ue = (Cdd/47r)[l — 3cos^((?£)]/^"^, multiphcd by the 
densities and Uj in the two sites, where 9£ is the angle between £ and the orientation 
of the dipoles. The ratio between the total on-site interaction U ^ Us + C/dip and the 
nearest neighbor dipolar interaction J/nn determines much of the physics of the system. 
It can be varied by tuning the on-site dipole-dipole interaction Udip from negative to 
positive by changing the vertical confinement, or by changing the s-wave scattering 
length via a Feshbach resonance, as recently demonstrated with Chromium atoms [70] . 

Due to the anisotropic character of dipole-dipole interaction, a much richer physics 
is to be expected with respect to the case of only repulsive long-range interactions. By 
changing the optical lattice strength, the transverse confinement, and the orientation 
of the dipoles, the parameters of the Bose-Hubbard Hamiltonian can be tuned over 
a wide range, the interaction can be made positive or negative and the tunneling 
parameter anisotropic. Exploiting all those degrees of freedom, one can scan in a 
single system checkerboard or striped ground states, and eventually collapse [269]. 

In [281], the case of Chromium, including the real spinor character of the 
atoms, has been considered. The dipole-dipole interaction has been treated 
as a perturbation on top of the dominant spin-dependent contact interaction. 
Using a mean-field treatment with a trial wavefunction beyond on-site product 
wavefunction, the quantum phases have been identified and in particular an 
antiferromagnetic-ferromagnetic first order transition occurring simultaneously with 



i All apart from [276]. 
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the MI-SF transition. 

The existence and observabihty of the above mentioned quantum phases require a 
relative strength of the long-range dipole-dipole interaction not too small compared to 
the zero-range one. This can be achieved by reducing of the s-wave scattering length, 
as demonstrated in [70]. However, the absolute energy scale has to be compared with 
recombination losses over the time-scale of the experiments and finite temperature 
effects. This might make stronger dipolar interactions desirable and polar molecules 
the optimal candidates for the realization of this kind of physics. Recently, as already 
mentioned in section 3.1, the difficulties of creating heteronuclear polar molecules in 
deeply bound vibrational states have been successfully overcome [24, 25, 26] and a 
high-space density gas of polar molecules in their ground vibrational state has been 
obtained [25, 282]. These achievements open a new era towards quantum degenerate 
molecular gases of strongly interacting dipoles. 

The search for supersolid and other exotic phases in cold atomic system in 
optical lattices does not however restrict to the case of dipolar atoms, Rydberg atoms 
or polar molecules [283, 284, 43]. Other optical lattice systems are relevant, like 
Bose-Bosc mixtures [285, 286, 287], Bose- Fermi mixtures [288, 287, 289], Fermi- Fermi 
mixtures [290], confined attractive [290, 291] and repulsive fermions [292], and bosonic 
gases in frustrated (triangular) lattices [293, 294, 295, 296], as well as extended Bose- 
Hubbard Hamiltonians designed from underlying contact-interaction systems using 
proper laser excitations involving higher bands [279]. 

9.5. Metastable states of dipolar lattice gases 

Beyond the richness of ground states discussed in the previous section, in the presence 
of long-range interactions, metastable insulating states are also predicted [297] . Similar 
physics appears also in the case of Bose-Bose mixtures, where local minima of the 
energy landscape indicate the presence of quantum emulsion states, i.e. metastable 
states characterized by microscopic phase separation, finite compressibility in absence 
of superfluidity, thus with macroscopic properties analogous to those of a Bose 
glass [298, 299, 300]. 

To determine the existence of the insulating states at fractional filling factor and 
the metastable states in a gas of dipoles in an optical lattice, one has to apply exactly 
the same criteria defining the Mott insulating states in the case of on-site interaction 
only. A crucial difference is that for non uniform atomic distributions in the lattice, the 
energy of particle- hole excitations is site-dependent. For nearest-neighbor interaction 
and zero tunneling, the checkerboard ordering of the atoms (see figure 31 (GS)) is 
the ground state in the range of chemical potential < /it < 4C/nn- In a similar 
way, an "elongated-checkerboard" ordering of the atoms (as shown in figure 31(1)) 
at J = is stable against particle-hole excitations in the range of chemical potential 
J7nn < M < 3/7nn, t^NN being the first nearest-neighbor interaction (and neglecting 
for the sake of simplicity all following ones). 

The phase diagram of the system for a 4 x 4 elementary cell and different cut- 
off of the interaction range is shown in figure 31. In the presence of long-range 
interactions, there are insulating lobes corresponding both to ground and metastable 
states, characterized by integer and fractional filling factors and a non uniform 
distribution of the atoms in the lattice (see the caption of figure 31 for details). 

The filling factors allowed and the metastable configurations clearly depend 
strongly on the cut-off range of interactions. Due to the strongly decreasing 
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Figure 31. (a,b,c) Phase diagram with a range of the dipole-dipole interaction 
cut at the first, second and fourth nearest neighbor respectively. The thick fine is 
the ground state and the other lobes correspond to the metastable states, the same 
color corresponding to the same filling factor. In (a,b,c) the ground state filling 
factors are multiples of 1/2, 1/4 and 1/8 respectively. In (a,b,c) the metastable 
state filling factors are m/4, m/8 and m/16 respectively, (Vm 1). Metastable 
configuration appearing at the first nearest neighbor (I), and second (Ila-IIb), 
and the corresponding ground state (GS); the metastable states remain stable for 
all larger ranges of the dipole-dipole interaction. The above phase diagrams are 
calculated for C//!7nn = 20. This value of the dipole-dipole interaction is much 
stronger than the one currently available in systems of Chromium atoms, where 
f^/f^NN ~ 400 (for e^d =~ 0.16 and spherical localisation of the bottom of the 
potential well at s = 20Eji, where Eji is the recoil energy at A = 500 nm). 

behaviour of the dipole-dipole interaction, in most theoretical approaches the 
interaction range is cut-ofF at few nearest-neighbors. This influences the results at 
very small particle or hole densities (in the particle-hole duality case of strong on-site 
interaction, as shown in figure 31), but only slightly changes the insulating part of the 
phase diagram at densities close to half filling [301]. This conclusion is confirmed by the 
results obtained in [302] for dipolar atoms in a ID optical lattice. Taking into account 
the infinite-range dipole-dipole interaction, one finds a Mott lobe in a given range 
on ji for each rational filling factor, but the dominant lobes are those corresponding 
to filling factors with smaller denominator. The influence of the extended range of 
interaction on the superfluid and supersolid part of the phase diagram deserves further 
investigation. 

This phase diagram is confirmed by the imaginary and real time evolution of the 
system. Depending on the initial conditions, the imaginary time evolution converges 
to a different metastable configuration. In the real time evolution, the stability of 
those metastable configurations is reflected into typical small oscillations around the 
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corresponding local minima of the energy landscape. 

The stability of the mctastablc states has been studied with a path integral 
formulation in imaginary time [303] , which can describe the tunneling below a potential 
barrier (instanton effect). This analysis suggests that the metastable configurations are 
very stable when many sites must invert their population to reach another metastable 
state. However, especially in larger lattices, two metastable configurations might 
differ just by the occupation of few lattice sites. This, and the corresponding small 
energy differences, should be carefully taken into account in a realistic analysis at 
finite temperature aimed at describing the experiments. 

Because of the presence of those very many metastable states, in an experiment 
it will be very hard to reach the ground state or a given metastable configuration. 
It was checked that the presence of defects is strongly reduced in the result of the 
imaginary time evolution, when a local potential energy following desired patterns is 
added to the optical lattice. One can use superlattices in order to prepare the atoms 
in configurations of preferential symmetry. This idea is presently pursued by several 
experimental groups. Note that the configurations obtained in such a way will also 
remain stable once the superlattice is removed, thanks to dipole-dipole interaction. 

Spatially modulated structures can be detected via the measurement of the spatial 
noise correlations function of the pictures produced after expansion [280, 256, 257]. 
Such a measurement is in principle able to recognize the defects in the density 
distribution, which could be exactly reconstructed starting from the patterns in the 
spatial noise correlation function. However, the signal to noise ratio required for single 
defect recognition is beyond present experimental possibilities, where the average over 
a finite number of different experimental runs producing the same spatial distribution 
of atoms in the lattice is required to obtained a good signal. 

In view of the possible application of such systems as quantum memories one 
should be able to transfer in a controlled way the system from one configuration 
to another [301]. The real time evolution of the system was studied with varying 
appropriately the Hamiltonian parameters (tunneling coefficients and local chemical 
potentials) and it was shown that it is impossible to map this problem onto a simple 
adiabatic transfer process. This is due to the fact that, in spite of the modification of 
the lattice parameters, the metastable states siuvivc unchanged till the point where 
the stability condition is not fulfilled anymore. The transition to any other state in 
then abrupt. A way around this problem is to push the system into the superfluid 
regicm and then drive it back into a different insulating state. The transfer between two 
metastable states turns out to be a quantum controlled process, where the Hamiltonian 
parameters must be controlled with very high precision to obtain the desired result. 

9.6. Bilayer and multilayer dipolar lattice gases. Interlayer effects 

In the previous subsections, we have discussed the case of a single 2D optical lattice 
layer. In present experiments, usually 2D geometries are created as a series of pancake 
traps by means of a very strong ID optical lattice in the perpendicular direction, which 
provides strong confinement and completely suppresses tunneling in that direction. In 
the presence of long-range interaction, in order to isolate each layer, one should also 
reduce the interaction between the different layers making the distance between the 
different layers much larger than the lattice spacing in the 2D plane. This can be 
achieved e.g. by creating a ID lattice in the perpendicular direction with two laser 
beams intersecting at a small angle 6, which increases the lattice spacing in the third 
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Figure 32. (a) Phase diagram for a system of two parallel ID lattices in 
the presence of onsite intra-wire interaction U and nearest neighbor inter-wire 
interaction W , for W/U = —0.75; white represents 2SF, gray MI, and black PSF 
(figure taken from [305]). (b) Effective MF phase diagram for a system of two 2D 
optical lattice layers in the presence of onsite intra-layer interaction U , nearest 
neighbor intra-layer interaction (7nn and nearest neighbor inter- wire interaction 
W , for W/U = —0.95 and U-mn/U = 0.05; the white regions inside the lobe are CB 
insulating with single or double site occupancy, the gray shaded region represent 
the PSS and the lower white region is PSF, as indicated by the arrows. The red 
line indicates the estimated limit of validity of the effective MF treatment [310]. 



direction to diD = (A/2)/ sin(6'/2). Alternatively, by using two different wavelengths 
for the ID and 2D lattices, one can make the ID distance larger or even smaller that 
the 2D lattice spacing. This might turn useful in cases where inter-layer interactions 
do not have to be suppressed, but on the contrary are exploited to generate novel 
effects. 

The case of two parallel ID optical lattices without tunneling among the two 
wires has been considered in [304, 305]. In that work, the polarization of the dipoles 
is chosen such that only onsite intra-layer interaction and nearest-neighbor attractive 
interaction between the layers exist. Such an attractive inter-layer interaction leads to 
the formation of a pair superfluid (PSF) [306, 307, 308, 309], i.e. a superfluid phase 
of pairs, composed by two atoms at the same axial position but in different wires. 
In the PSF phase, only simultaneous hopping of atoms in the two wires is involved. 
Due to the direct JVIott to PSF transition the lowest excitations of the JVIott state 
are not, as usually, given by particle-hole excitations, but rather by the creation and 
destruction of pairs. This change in the nature of the Mott excitations leads to a 
significant deformation of the Mott-insulator lobes, and may even induce a re-entrant 
shape of the lobe at small hopping. 

The same scenario is expected for two 2D optical lattice layers, based on the 
mapping of this problem to the one of bosonic mixtures in 2D lattices [306, 307, 308, 
309]. However, the physics of dipolar atoms in layered 2D optical lattices is even richer, 
because of the long-range intra- and inter-layer anisotropic interactions. The case of 
dipoles pointing perpendicular to the 2D lattice plane, generating repulsive nearest- 
neighbor intra-layer and attractive nearest-neighbor inter-layer interactions, has been 
recently investigated [310]. A mean-field treatment of the effective pair Hamiltonian 
provides clear evidence of the existence of a pair-supersolid phase (PSS), which arises 
from two-particle and two-hole excitations on top of the checkerboard-like insulating 
phase at half-integer filling factor (see figure 32, right panel). 

When tunneling between the two layers is not completely suppressed (in the 
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specific case, two uniform 2D layers without lattice), beyond the superfluid and pair- 
supcrfluid phases, a phase transition towards a maximally entangled state, where all 
particles populate either one layer or the other, has been shown [311]. 

In [312] the full 3D geometry for dipoles pointing along the perpendicular 
z direction has been considered. For inter-layer attraction and positive on-site 
interaction, layered phases where the density distribution is the same on all lattice 
layers exist. For negative on-site interaction, instead, modulation of the density with 
period 3 along z has been found, reminiscent of the structure of high-Tc cuprate 
superconductor. In this work there is no evidence of the pair superfluid phase, because 
of the mean-field approach used. 

9. 7. Tailoring interactions with polar molecules 

Thanks to their promisingly large dipole moments (typically of the order of a few 
Debye), polar molecules have been the subject of extensive theoretical investigation 
centered in the possibility of tailoring the shape of their interaction [313]. 

Polar molecules prepared in a mixture of rotational states interact through long- 
range dipole-dipole interaction even in the absence of an external electric field. The 
possible Mott phases present different ordering, depending on the preparation of 
the initial superposition. When the Mott state is melted the superfluid state can 
interpolate between homogeneous and antiferromagnetic ordering or phase separate 
depending on the Hamiltonian parameters [284]. 

On the other hand, appropriate static and/or microwave fields can be applied to 
design effective potentials between two molecules in their electronic and vibrational 
state [43, 314]. The mechanisms to induce tailored interactions relies on the simple 
rigid rotor Hamiltonian (see Appendix B), providing the low energy rotational states 
of the molecules. Such rotational states can be coupled by static or microwave 
fields to design long-range interaction potentials, in contrast to the typical van der 
Waals interaction in absence of external fields. In 2D, a static electric field induces a 
typical repulsive dipolar potential scaling as l/r^, which leads to crystallization. An 
additional microwave field allows further shaping of the potential, even allowing an 
attractive part. As previously discussed, in the presence of such induced interactions, 
the system can be driven through crystalline, superfluid and normal phases [43, 314]. 

Further investigation has shown that the static electric field and the microwave 
field dressing the molecular rotational levels can be chosen in such a way to obtain 
dominant three-body interactions [315]. Hamiltonians with many-body interactions 
have been studied in the contest of non abelian topological phases (like Pfaffian 
wavefunction accounting for the quantum Hall effect or systems with a low energy 
degeneracy characterized by string nets), multiple species in frustrated lattice 
topologies, ring exchange models (like the one responsible of the nuclear magnetism in 
Helium 3), or undoped high-Tc compounds and cuprate ladders. In typical condensed 
matter systems many-body interactions are rarely dominant and polar molecules 
provide a setting where they can be controlled and designed independently from two- 
body interactions. 

A Hubbard model including an unconventional three-body interaction term 
T^i^j-tk WijkniTijnk can be readily obtained. The Hamiltonian parameters for two- 
and three-body interactions depend explicitly on the applied fields. The three-body 
interaction is intrinsically long-range due to the underlying dipole-dipole interaction. 

The 2D and ID quantum phase diagrams for strong three-body interactions have 
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Figure 33. (a) Square lattice in 2D witli nearest-neighbor orientation-dependent 
Ising interaction along x and z. (b) Two staggered triangular lattices with nearest- 
neighbor oriented along orthogonal triads. Figure taken from [283]. 

been recently investigated. In 2D [316], a rich variety of solid, supersolid, superfluid 
or phase separated phases are encountered. The several solid phases at fractional 
filling factor evolve, upon doping, into corresponding supersolid phases with complex 
spatial structures. In particular, the checkerboard supersolid at filling factor 1/2, 
which is unstable for hardcore bosons with nearest neighbor two-body interaction, is 
found to be stable in a wide range of tunneling parameter. In ID [317], quantum 
Monte Carlo simulations have shown that strong three-body interactions give rise to 
an incompressible phase at filling factor 2/3, which presents both charge density wave 
(CDW) and bond (BOW) orders. At the same time, they have ruled out the solid 
phases at filling factors 1/2 and 1/3 predicted by Luttinger theory [315]. The solid 
phases at filling factor 1/2 is found only in the presence of additional nearest and 
next-nearest two-body interactions. These can be either CDW or BOW depending on 
the intensity of the two-body corrections. Instead, at filling factor 1/3, the system is 
always superfluid. 

By tuning the optical potential parameters and by means of external electric 
and magnetic fields, one can induce and control the interaction between spin states 
of neutral atoms in an optical lattice and engineer quantum spin Hamiltonians [318, 
319, 17]. Such proposals are aimed to the study of a variety of quantum phases, 
including the Haldane phase, critical phases, quantum dimers, and models with 
complex topological order supporting exotic anyonic excitations. Including a spin 
degree of freedom in addition to the rotational degrees of freedom of polar molecules, 
spin models for half-integer and integer spins with larger coupling constants can be 
obtained [283, 320]. The main ingredient of these proposals is the dipole-dipole 
interaction: it couples strongly the rotational motion of the molecules, it can be 
designed by means of microwave fields (as explained above), and it can be made spin- 
dependent, exploiting the spin-rotation splitting of the molecular rotational levels. The 
final goal is to reproduce models with emergent topological order, robust to arbitrary 
perturbations of the underlying Hamiltonian, and hence suitable for error-resistant 
qubit encoding and for quantum memories. 

For spin 1/2, it has been explicitly demonstrated how to construct two highly 
anisotropic spin models [283]. The first model (see figure 33(a)) is a 2D spin- model 
with nearest-neighbor orientation-dependent Ising interactions. It has ideally a gapped 
two-fold degenerate ground subspace with zero local magnetization, which guarantees 
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immunity to local noise [321]. The second model (see figure 33(b)) takes place on 

two staggered triangular lattices, equivalent to Kitaev's honeycomb model [322]. In 
appropriate regimes, this model provides a ground state which encodes a topologically 
protected quantum memory. 

9.8. Self- assembled structures 

The long-range character of the dipole-dipole interaction allows the formation of self- 
assembled structures. Different situations have been the object of recent studies, 
ranging from the formation of chains of polar molecules in ID optical lattices driven 
by the attractive part of the dipole potential, to the appearance of crystal ordering in 
ID or 2D systems driven instead by the repulsive character of the interaction. In this 
section, we summarize the main results and possible applications. 

The case studied in [323] considered a stack of 2D layers created by a strong 
ID lattice. The dipoles are pointing perpendicularly to the 2D layer, such that the 
interaction is repulsive in each layer and attractive for atoms located on top of each 
other in different layers. Collapse in the perpendicular direction is prevented by the 
complete suppression of tunneling by the strong ID lattice. It has been found that the 
attraction in the perpendicular direction is responsible for the formation of chains of 
dipoles, where the longest chain are energetically favorable, while the shortest chain 
are favored by entropy. When temperature is decreased, condensation in the longest 
chain takes plac;e. It has been pointed out that the physics is similar to the physics of 
rheological electro- and magneto-fluids. 

In the 2D homogeneous system for dipoles pointing perpendicular to the 
plane, the repulsive interaction is responsible for the formation of hexagonal crystal 
ordering [324, 43, 325, 326], similar to the formation of Wigner crystals for electron 
or ions interacting via the Coulomb repulsion [327, 328]. However two important 
differences are found, namely that the transition to the crystal phase in dipolar system 
happens at high density (instead that at low density as for Coulomb systems), and 
the spectrum shows two phononic linear branches with different slopes (rather than 
the dispersion u) ~ g^/^ typical of Coulomb crystals). 

The formation of the crystal ordering has been identified as a first order phase 
transition [43, 326] which appears with a delta-peak at the inverse lattice spacing in 
the structure factor. The transition happens for = 37 ± 1 [326], where r,, stands 
for the ratio between interaction energy and kinetic energyff . Contrary to Coulomb 
systems, increases with density, that is why Bose-Einstein condensation is found 
at low density and the dipole-dipole dominated crystalline phase at high densities. 
The quantum melting transition should be in reach for the typical parameters of polar 
molecules. 

In [325] the excitation spectrum has been shown to develop a roton minimum 
as the density increases. However, the question about the possibility of having a 
supersolid in such systems is still open due to difficulties in getting reliable measures 
of the superfluid density [326]. In [114] a greater stability of supersolid phases has been 
conjectured for multi-layer systems, since the roton softening occurs at wavelengths 
much larger than the layer width, which should prevent global collapse. 

In [329], it has been proposed to exploit the protection against short range 
collisions provided by the crystalline ordering to use those self-assembled dipolar 
crystals as high-fidelity quantum memories. Furthermore, the self-assembled crystal 

ft This result is consistent with the other results of ra = 32 ± 7 [43] and Ts = 30 [325]. 
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could be used as a floating lattice structure for a second species of atoms or 
molecules [330]. The advantages of this proposal are that the lattice spacing can 
the be tuned down to the hundred nanometers scale and that the second species is 
subject both to its own Bose-Hubbard coupling and to the coupling to lattice phonons. 

In a system of two parallel 2D layers at distance I and dipoles pointing 
perpendicular to the layers, one finds intra-layer repulsion and inter-layer attraction 
(repulsion) for dipoles separated by a distance r < \/3Z (r > \/3Z). It has been 
shown [331] that the hexagonal crystalline order is preserved in both layers, and as a 
function of the layer distance /, one goes from two independent crystals to a crystal 
of paired dipoles. Moreover, since the melting temperature is not monotonic in I, a 
solid-liquid-solid transition takes place for increasing I at fixed temperature. 

Furthermore, the formation of ordered patterns in ID systems has been recently 
investigated in [332, 333, 334, 335, 34]. In [333, 334], it has been pointed out that 
one-dimensional dipolar gases present strongly correlated phases beyond the strongly 
correlated Tonks-Girardeau regime. The crossover from the superfiuid state to the 
ordered state takes place for increasing densities (like in 2D) and appears in the 
structure factor as additional peaks at the inv(;rse lattice spacing. In the whole 
crossover, the gas preserves a Luttinger-liquid behaviour, since no roton minimum 
nor long-range order are found. 

In [34] the stability of such ordered structures with respect to the transverse 
confinement has been investigated, in analogy to what is known from Coulomb Wigner 
crystals [336]. By weakening the transverse confinement, or equivalently increasing 
the density or the strength of the dipolar interaction, a smooth crossover to a zigzag 
chain and to structures formed by multiple chains is predicted. Quantum fluctuations 
smoothen out these transition, which are respectively first and second order phase 
transitions in a classical model, and also completely melt the crystal for low values of 
density or dipolar interaction. 

10. Outlook 

An epilogue, in the disguise of wrapping up the past, 
is really a way of warning us about the future. 
J. W. Irving, The World According to Garp. 

10.1. From Chromium to heteronuclear molecules to Rydberg atoms 

The experimental realization of dipolar gases was first obtained with atomic 
Chromium. Despite the large spin magnetism in Chromium, the dipole moment is 
still small, classifying the dipolar interaction as a weak one, where the length scale 
of the interaction is much smaller than the interparticle spacing. In this regime we 
will probably see in the next future experiments on self-organized collective structures 
close to the instability of the gas. In spinor gases first experimental steps have been 
done in this direction [208], as shown in section 8. 

Rapid progress has been made in the creation of heteronuclear molecules (see 
section 3.1) in their vibrational ground state [25, 26]. These experiments still have 
to go some way to degeneracy; however given the speed of the development we can 
expect this step to be taken in the next future. These gases — when the molecules 
are prepared in the rotational ground state — then cover the whole range of dipolar 
interaction strengths — from weak to strong interactions — as their electric dipole 
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moment can be tuned via a external electric DC field. Even the three body interaction 

can be controlled independently of the two body interaction by the use of microwave 
fields [315]. These systems will therefore provide a rich toolbox for quantum simulation 
of spin systems. 

Even dramatically larger dipolar interactions between Rydbcrg atoms have 
also become available in the regime of quantum degenerate gases since the first 
Rydberg excitation of Bose-Einstein condensates [49]. As the energy spectrum of 
Rydbcrg atoms involving a quantum defect is very similar to the level structure of a 
heteronuclear molecule, many of the proposed techniques to mix rotational states to 
tailor the dipolar interaction can be directly applied to Rydberg atoms. Despite their 
limited lifetime the huge size of the dipolc moments typically 1000 times larger 
than for heteronuclear molecules — opens up a whole new class of experiments on 
long range interacting spin systems. Recently the mapping of a long range interacting 
spin system onto a frozen Rydbcrg gas was realized successfully and universal scaling 
of that strongly interacting system could be determined experimentally [337]. 

In terms of particle numbers this experiment was done in a regime where exact 
quantum calculations have no access. Therefore the quantum simulation of large and 
strongly interacting spin systems by mapping to other equivalent quantum systems 
— just like Richard Feynman was proposing — is starting to become reality. Dipolar 
gases will enrich this effort substantially. 

10.2. Dipolar gases and trapped ions 

One of the most promising and fascinating ideas that will allow to investigate quantum 
systems baring a very close similarity to ultracold atomic gases comes from somewhat 
unexpected directions. Namely, in the recent years, there has been an enormous 
interest in ideas that could lead to use trapped ions for simulating quantum many 
body systems. 

Ch. Wunderlich had in 2001 [338, 339] the idea of using inhomogeneous magnetic 
fields and long wave-length radiation to solve the problem of individual ion addressing 
in ion traps one of the most important obstacles toward the implementation of 
quantum information processing with ultracold ions. As a byproduct Mintert and 
Wunderlich obtained the result that internal degrees of freedom of ions behave as 
coupled pseudo-spins, where the coupling is mediated by the phonons in the trap. This 
idea has been fully put forward by Porras and Cirac [340, 341], who have argued that 
trapped ions can be as attractive systems for quantum simulations as ultracold atoms 
or molecules in optical lattices. Effective quantum spin systems that one can achieve 
with ions range from linear chains to 2D self-assembled, or optically prepared lattices, 
with very precise control of the parameters. The most interesting aspect in the present 
context is that the interactions mediated by phonons are of long range character, and 
in fact typically decay with distance as 1/r^. The typical spin Hamiltonian one can 
simulate with ions has the form 

rr _ ■'■ \ ^ rz qz qz ^ \ ^ rx ax qx ^ \ ' jV qD qH 
ij ij hj 

-J^iHfSf + HfSy + HfS!), (10.1) 

i 

where the overall sign and magnitude of J^j^'^' can be controlled, whereas the spatial 
dependence is J^j^'^ oc l/\ri — rj\^. When = J" these models correspond nearly 
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exactly to hard core Bose gases with dipole-dipole interactions, and in the presence of 
transverse and longitudinal "magnetic" fields. The only difference is that the tunneling 
now has a non-local character, and its amplitude also decays as slowly as 1/r^, since 
the phonons allow for such long-range interactions. 

The ingenious idea of Porras and Cirac has been followed in many directions. 
Phonons themselves may exhibit fascinating many body quantum physics and undergo 
Mott insulator-superfluid transition [342]. One can realize mesoscopic spin-boson 
models with trapped ions [343] , or quantum neural networks models and fault resistant 
quantum computing [344, 345] . Recently, the idea of combining lattices (either optical 
ones consisting of microtraps [346], or lattices employing surface electrodes [347]) 
with ion traps has led to fascinating proposals for realization of antifcrromagnetic 
spin models in a triangular lattice with adjustable couplings, where the different Neel 
orders can be realized [346] . As the parameters interpolate between these Neel orders, 
the system is expected to enter into quantum spin liquid states of various types. These 
predictions have been very nicely supported in [346] , where exact diagonalization was 
compared with one of the first nontrivial applications of the so-called Pair Entangled 
Projected States (PEPS) method. 

All these proposals are not inventions of theoreticians: many of the leading 
experimental ion trappers groups are working on quantum simulations with ions, and 
the first spectacular experimental results have been obtained already [348]. 

10.3. "Dipolar art" 

To conclude this review, we note that the excitement about the advances in dipolar 
gases is not only widespread among the experts in the cold atom and condensed matter 
community: to our amazement, it has also inspired recently the artist Brigitte Simon 
to work on "dipolar art". This expert in artwork made of glass has for example 
designed many church stained glass windows. After discussing the physics behind the 
images of a collapsing dipolar BEC with some of us, she started to work on a glass 
window based on our measurements. Let us quote her: "I was fascinated when I 
discovered an article in the Frankfurter AUgemeine Zeitung in August 2008 [about the 
dipolar collapse, reproducing pictures from figure 19] . Reading it over and over I tried 
to understand what a miracle had been created at Stuttgart University. Marveling at 
the photos I felt that there is Art in Science, and this is what I tried to capture" . 
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Figure 34. A stained glass window, composed of 524 one-inch squares, each 
selected from a wide variety of colors and precisely cut out of mouth-blown 
Genuine Antique Glass from Bavaria, reproducing the shape of a collapsed dipolar 
BEG (compare with figure 19). The squares are framed with copper foil and 
soldered front and back. Figure courtesy of B. Simon. 
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Appendix A. Fourier transform of the dipolar interaction 

In this appendix, we sketch the main steps of the calculation of the Fourier transform 
(2.5) of the dipole-dipole interaction. Using spherical coordinates {r,d,ip), with the 
polar axis along k and the dipole moment in the y = plane (making the angle a 
with fe), one has 

U^,^k)^^JJJ e-^^ - ' 1 - 3(sina sing cos y. + cosa cos 0)^ ^.^ ^ ^^^^^ 

After integration over tp and the change of variable x = cos 9, we obtain 

C^(fc) = 9^ ^ y ' e-*'="^(3 cos^ a - 1) (1 - 'ix') dx (A.2) 

where 6 is a cutoff at small distance introduced to avoid divergences at this stage. The 
integration on x is straightforward, and gives 

„ „ 2 X f°° f sinu 3cosu 3sinu\ , , . „^ 
f/dd(fc) = Cdd(l-3cos2a) / (^r + 3 4— du,(A.3) 



J/cb \ " " " / 

with u — kr. The last integral can be calculated by parts and has the value 
[A:6cos(fc6) — sin(fe6)]/(A;6)'^. We can now let the cutoff b go to zero; the last integral 
then approaches —1/3 and thus we finally get the expression (2.5) for the Fourier 
transform of the dipole-dipole interaction. 
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(a) (b) (c) 
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Figure Bl. (a) Rotational spectrum of a diatomic molecule in zero field, (b) 
Dependence of the first energy levels on the applied electric field £. (c) the 
ground state average dipole moment (dz) in the laboratory frame as a function of 
the applied field £. 



Appendix B. Stark effect of the rigid rotor 

In this appendix, we briefly recall basic results (see e.g. [314]) concerning the behaviour 
of a spinless diatomic molecule, modeled as a spherical rigid rotor, in an electric field 
£, with an emphasis on the dependance of the average electric dipole moment (in 
the laboratory frame) on the applied field. We assume that the molecule is in its 
electronic and vibrational ground state, and that the electronic ground state is a "'^S 
state (like e.g. in the case of bi-alkali molecules). For the sake of simplicity, we also 
neglect the hyperfine structure, although this is an issue relevant to experiments. The 
hamiltonian for a rigid rotor reads 

Hrot = Bj\ (B.l) 

where J is the molecule angular momentum operator (in units of h) and B the 
rotational constant, linked to the equilibrium internuclear distance r and the reduced 
mass mj- by the relationship B — h"^ / {2m^r'^); its typical order of magnitude is 
B/h ~ 10 GHz. The eigenstates of (B.l) are the angular momentum eigenstates 
|J, TOj) with energy BJ{J + 1), and are 2 J + 1 times degenerate. Figure Bl(a) 
represents the first few eigenstates of (B.l), with energies 0, 25, 65. . . 

This molecule is supposed to have a permanent dipole moment d in the molecular 
frame. Then, in the presence of an external field E — Ebz, the Stark hamiltonian of 
the molecule reads: 

H = H,ot -d E = H,ot - dS cos 0, (B.2) 

where 9 is the angle between z and the internuclear axis. Figure Bl(b) represents 
the first eigenstates of the hamiltonian (B.2), diagonalized numerically, as a function 
of £. The interaction with the electric field lifts the degeneracy between levels 
having different values of \mj\. From this Stark map, the average dipole moment 
{dz) = d{cos6) for the ground state {(jjo) is obtained (via the Hellmann-Feynman 
theorem) as 

{dz) = - {(jJo 



dH 



d£ 



' d£ 



(B.3) 



CONTENTS 



64 



where Eq{£) is the ground state energy. The dipole moment is plotted on figure Bl(c). 
One observes that {dz) increases linearly at small £ (more precisely, one has 
{dz)/d ~ d£/{3B) for small fields), and tends asymptotically for d£ B towards 
its saturated value d, although relatively slowly, as one can show that to leading 
order, (4)/o! - 1 - y/B/{2d£) for d£ :$> B [349]. For typical values d - 1 D and 
B/h 10 GHz, the electric field strength corresponding to d£ ^ B is on the order 
of 10^ V/cm, which, from the experimental point of view, is accessible in a relatively 
easy way. 
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